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1. INTRODUCTION 
THE well-known equations due to Blasius 


5 iG w'dz = Fy — iFy (1) 


~§ § w'? zdz=M+iN (2) 


give the force and moment exerted by a perfect, incompressible, irrotational 


fluid upon material bodies immersed in it, in the case of steady motion in 
two dimensions. The real part of (2) gives the moment, and the imaginary 
part, the quantity N which is called the virial of the forces plays an im- 
portant role when we study the changes that occur when the forces are rotated 
by a certain amount about their points of application without change of 
magnitude. For such a change the new resultant has a line of action also 
turned by the same angle but without change of magnitude. As the angle 
of rotation is changed from 0 to 27 we will get a family of straight lines 
having a certain curve as its envelope, and the point where a line of action 
touches the envelope or the point at which the resultant force rotates for an 
infinitesimal rotation of all the components is called the Hamilton Centre. 


In particular, for a system of forces F, (X,, Y,) acting at (x,, y,) in the plane, 
the moment about the origin is M= Y (x, Y,—y,X,) and N=J (x,X,+y,Y,) 


is the virial referred to the origin. It is well known that the Hamilton 
Centre is the intersection of the resultant 


M = xY — yX (3) 
where X = DX,, Y =ZY,, and the line perpendicular to (3) given by 
N = xX + yY. 4 


If the Hamilton Centre coincides with the origin we get M = 0, and N = 0. 
In other words, the Hamilton Centre C is the point for which both the moment 
and virial vanish. 
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The point C has a significance similar to that of the metacentre P. This 
point is the point of intersection of the original resultant, and the new line 
of action obtained when we consider not only a change in the direction of 
the pressure but also the consequent simultaneous change of magnitude in 
accordance with the hydrodynamical equations. The point C is obtained 
when the magnitude of the pressure force at any element does not change. 
It is well known that for stability the centre of gravity G should lie below P. 
In the case of a real disturbance under flying conditions, it may be that at the 
very moment a small angular disturbance occurs the pressure at a point does 
not change in any way. In such a case we should consider the. point C 


instead of P. Thus, in general, one has to ensure for stability that G must 
lie below both P and C. 


The general theory sketched above is due to von Mises (1941, pp. 43-47), 
but nowhere in the literature does one find problems relating to the 
Hamilton Centre dealt with in any specific case. Thus in the case of the 
simple wing profile, i.e., a cylindrical wing of infinite span with a cross- 
section containing one singular point, expressions for the lift force and mo- 
ment have been derived, and we have the well-known theorem of von Mises 
that the first axis is the directrix of the metacentric parabola. But no results 
relating to the virial have however been obtained. 


It is the purpose of the present paper to derive some theorems relating 
to the virial; and the Hamilton Centre, in the case of the simple wing profile 
and particular cases of this like the Joukowski and Karman-Trefftz profiles, 
and the S-profile due von Mises. We also derive corresponding results 
relating to stability in these special cases. A parabola termed the virial 
parabola is introduced, and a theorem analogous to that of von Mises has 
been obtained. Further a relationship is established between the Hamilton 
Centre and the centroid of circulation (Milne-Thomson, 1952, pp. 120-22). 


2. REVIEW OF RESULTS IN THE GENERAL CASE 


Considering the effect of uniform flow specified by the velocity vector 


Vv (or — Ve-*) at z= co upon immersed bodies, let W (z) be the complex 
potential of the two-dimensional (vortex free) motion of the liquid of 
density p. W’ (z) is given by its Laurent expansion 


Wi) =A+54+ SH... (5) 


where it can be shown that A, B, C are linear in V, and that A = — Ve", 
B = I'/2zi, where I is the circulation. 
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Writing C=V (ce + ce*); B=Vi (—), cos a+b, sin a), 
2a = | V(b,2+ by?) |; cos B = — b,/2a, sin B = b,/2a, we get 
= — 2a Visin(a+ §). Using Blasius’ formula (1), we get the lift force 
L = pI'V = 4mapV? sin (a + 8) = Ly sin(a + 8) (6) 
in a direction normal to the direction of flow. Also frdm (2), taking the 
moment and virial relative to a point Zz), we have 


My, + iNy = -§ G (2 — %) w'?2 dz 


ee 5 (B? + 2AC) 2ai + 5 70 . 2AB . 2m (7) 


Using the properties of A, B, and C, it can be shown from this that if z» be 


taken to be the point F given by z, = — <! -18 the moment Mg, is independ- 


ent of the angle of attack. The point F is called the focus, and referred to 
F, we get 


My = — 2mpV?2 Gm (cz + c,e-*8) (8) 
Introducing hy = M,/ No (9) 
one can draw a line 1, in the 8-direction (called the direction of the first axis) 
at a distance hy from the focus F. If a parallel to V through this point 


> 

meets 1, at S, it can be easily shown that the normal to V through S gives 
the direction of the lift force L and the correct magnitude follows from (8). 
For varying a, S moves along 1, and L envelopes a parabola called the 
metacentric parabola, and the point of contact of L is the metacentre P. In 
the above figure where A, > 0, the metacentric parabola opens upwards, 
and if h, <0, it opens downwards. Thus the condition for stability, viz., the 
centre of gravity G lying above P is easier to reach in the case hy > 0 than 
in the case hy <0. 

We can find from (7) the value of N, and show that if we choose the 
point Gy given by z) = — - eB + ge’®, then Ny is independent of a. 
The point Gy may be called the virial focus, and referred to it we get 

No’ = 2npV? R (cz + cye~2**) (10) 
Introducing ky = Ny’ /Ly, (11) 
we plot the straight line m) perpendicular to 1) at a distance ky from Gp. 
Then as in the case of the moment, it can be shown that we get the line 


er 
Ny = 0 by drawing a normal to V through G, to meet my at T and taking 
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the line through T parallel to V. The intersection of M, = 0, and N, =0 
gives the Hamilton Centre C. 


3. RESULTS IN THE CASE OF THE SIMPLE WING PROFILE 


In the case of such a profile with one singular point B, the problem 
amounts to finding a transformation z’ = f(z) transforming the profile into 
a circle of radius a and centre M, leaving unchanged the region at infinity, 
and making the point B’ corresponding to B a stagnation point. If w(z) 


be the complex potential satisfying the required conditions, it can be shown 
that 


w! (2) = — Velie a , : + (Vat eis — Vei# k) - ++ (12) 


Comparing this with (5), we have A = — Ve"*, and B = I'/2zi as obtained 
previously. Also we have 


Cc, = —k, Cy = a. (13) 


The first axis is in this case the line B’M, i.e., the line joining the mapping 
point of B and the centre of the profile. Putting k in (13) equal to ce’, 
the inverse transformation z= g(z’) for large | z’ | transforms the circle 
| z' |= R into an ellipse whose major axis makes an angle y with the real 
axis, and this major axis is called the second axis of the profile. The lift force 
is again given by (6), and to find My we note that the point z) is now 
Bo 
= — e(27-8), This shows that that the second axis bisects the angle 
between the first axis and the line MF joining M to the focus F (z,), and also 
gives a construction to find F if we note that MF = c?/a. 
M, given by (8) now becomes My = 27pV?R ( —ike-**8), and since k = ce”, 
we have 
M, = 27pV2c? sin 2 (y — B) (14) 
and 
— SS A oe 
hy = .. sin 2(y — 8) (15) 
From (15) easily follows the theorem of von Mises that the first axis 
is the directrix of the metacentric parabola. Finally we can find the moment 
about M knowing it about F, and we get 


M = 2zpV*c? sin 2 (a + y) (16) 


It might also be noticed that axis II is a tangent to the metacentric 
parabola for, if FRJ be drawn perpendicular from the focus F to axis Il 





















Virial Problems Related to Simple Wing Profiles 57 


0 meeting it in R, and the first axis in J, we then have FR = RJ since axis II 
bisects the angle between FM and MJ. Hence R lies on the tangent at the 
yertex to the metacentric parabola, and this shows that axis II is a tangent 
to the parabola (Fig. 3). 
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. Fic. 1. Méetacentric Parabola, and lift force (4) > 0) 
le 
4. THEOREMS RELATING TO THE VIRIAL 
. We now come to the consideration of the virial in the case of the simple 
wing profile. Analogous to the metacentric parabola we can think of the 
", envelope of Ny = 0 as the virial parabola, and the point of contact, Q, of 
N, = 0 with this envelope as the virial centre (Fig. 2). Using (10) and 
4) putting here c, = — k = — c*e*” and c, = a® we get from (11) that 
e. ¢ 
‘ ky = 5 — 9, 0082 (y — B) (17) 
Let us find the distance of the directrix of the virial parabola from the centre 
is of the profile. This is obviously equal to FB + the resolved part of MF 
nt along axis I, i.e., equal to a + 2k) + c?/a cos 2 (y — 8), and using (17) this is 
equal to 2a. Hence we have 
6) THEOREM 1.—The directrix of the virial parabola is at a distance from 
ic the centre of the profile equal to the diameter of the profile. 
I If P and Q be corresponding positions of the metacentre and the virial 


centre, the tangents to the two parabolas at P and Q are invariably at right 
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angles intersecting at the Hamilton Centre C. Also from the properties 
of the parabola (Fig. 4), triangle FSP is similar to FOS, and the triangle 


Fic. 2. Construction for Hamilton Centre (4) > 0, Ky > 0) 


G,TQ is similar to Gj)KT. Hence FPS = OSF = a+ 8, and G,QT= 
G,TK =a-+ 8. Thus FP and G,Q make equal angles with SP and QT 


which are perpendicular directions. Hence FP and G,Q are at right angles, 
and: we have 


THEOREM 2.—The lines joining the metacentric focus to the metacentre, 
and the virial focus to the virial centre are at right angles. 


Next, we shall consider the centroid of the circulation given by 


fe ees (a Bp {a? et — c2e@i epee (18) 


(Milne-Thomson, p. 121). Let K, be the point diametrically opposite to 
B on the circle which transforms into the profile, i.e., K, is a point on the 
first axis such that MK, = a(Fig. 4). It is known that the locus of the 
centroid of circulation is the perpendicular bisector of FK,, and is called 
the third axis of the profile. Also the line of action of the aerodynamic force 
passes through the centroid of circulation. This point is therefore the inter- 
section of axis III with the tangent to the metacentric parabola at P, and 














Virial P oblems Related to Simple Wing Profiles 





Fic. 3. Axes I and II, and the Focus (hi, < 0) 


is denoted by C’ in Fig. 4. Since K, is a point on axis I which is the direc- 
trix of the metacentric parabola whose focus is F, the perpendicular bisector 
must be a tangent to the parabola. Thus the third axis is a tangent to the 
metacentric parabola, and the point of contact is denoted by C” in Fig. 4. 


The points S, C’, P and C all lie on the tangent at P to the metacentric 
parabola, and we will now find the distance CC’. Taking O for origin 
(Fig. 4), 1, for x-axis, and the axis of the parabola for y-axis, the equation 
to that parabola becomes 

x? = 4hoy (19) 


The tangent at P is perpendicular to the wind direction, and its slope is 
therefore equal to cot (a + 8). The angle a + 8 which is the angle between 
the direction of the wind and the first axis is called the absolute incidence, 
and we will denote it by fp». Putting cot 8) = ft, the point P becomes 
(2hot, Apt”), and the equation to the tangent at P is 
xt— y=hft? (20) 
It is easy to show that Q is the point (a + ky — kol®, hy + kot), and the 
tangent at Qto the virial parabola has the equation 
x+ yt =at+ky t+ hot + kot? (21) 
The intersection of (20), and (21) gives the Hamilton Centre C as 
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2 
We see from Fig. 4, that AK, = MK, — AM =a— A cos 2 (y— ) 


= 2k, from (17). Hence K, is the point (2k), — ho) and axis III which is 


the perpendicular bisector of FK, has the simple equation 
=] 


9 
xX 
zt 





Fic. 4. Showing the Axes 
xky — ho = ky? (23) 
It is easily seen that this is a tangent to the metacentric parabola the point 
of contact being (2k, k»?/ho). The centroid of the circulation is given as 
the intersection of (20) and (23), and has the co-ordinates 


C’ (ky + Mot, Kot) (24) 
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From (22), and (24) we calculate CC’, and we have 


CC’ = a/V/(1 + #) = asin By (25) 
leading to 


THEOREM 3.—TZhe distance of the centroid of circulation from the 
Hamilton Centre is equal to the perpendicular from the stagnation 
point on the direction of the wind through the centre of the profile. 


This theorem gives an alternative construction for the Hamilton Centre. 
We mark MK, on the first axis such that MK, = a, then draw the third axis 
to meet M, = 0 at C’, and then make C’C = equal to the perpendicular 
from K, on the direction of the wind through M. 


By drawing suitable figures we can find the relative positions of C and 
P for the cases hg 20, kj) =O. Summarising these cases we have 


THEOREM 4.—If hy> 0 (or 8 < y), the Hamilton Centre is above or below 
the metacentre according as ky =0, and 


THEOREM 5.—if hy <0 (or B> y), the Hamilton Centre is above or below 
the metacentre according as k,= 0, provided in the latter case 
ky > asin®Bo, otherwise C lies above P in the case ky> 0 also. 


These theorems enable us to make statements about stability based on the 
relative positions of G and one of the points C or P. 


Eliminating ¢ in the co-ordinates of C given by (22) we can derive the 
locus of C, and since both the co-ordinates are rational functions of a para- 
meter, the locus will be an unicursal curve, a cubic as can easily be shown. 


The distance between the metacentre and the Hamilton Centre can be 
found from (22) and using the fact that F is (2hot, hot*), and we get 


CP = asin By + ky cosec By — My cosec By cot By (26) 
as can also be deduced geometrically by taking (see Fig. 4) 
CP = G,T + G,D — SP. 


Considering the question as to when the Hamilton centre lies respec- 
tively on one of the three axes of the profile, it is obvious that since axes II 
and III are both tangents to the metacentric parabola, the tangents will also 
be the directions of the lift forces for suitable a. Thus C will lie on these 
axes when the wind directions are perpendicular to them. As regards 
axis I the result is trivial since when the lift force is along the tangent at the 
vertex of the metacentric parabola, C will lie on axis I at infinity, 
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Y 








Fic. 5. S-Profile-Three Singularities 
5. SPECIAL PROFILES 


We will now consider the special profiles like the Joukowski, Karman- 
Trefftz and S-form profiles. For the Joukowski profiles we have y =0 


and ky = = (a? — c* cos 28) and if the camber be c/a, we have c <a, <a 


i.e., a >c and therefore a > c cos 8, i.e:, ky > 0. Also since y =0 we have 
B>vy, i.e., 4g <0. Thus, in the general Joukowski profile 4, <0 and 
k, > 0 and it can be shown that C lies above P. Although in this case the 
metacentric parabola opens downwards, and it is harder to make P lie above 


G for all angles of attack, the situation is more favourable with regard to 
C than with regard to P. 


The results are, of course, true for the particular cases of the symmetric 
Joukowski profile (8 = 0), and the still more particular case of the straight 
slit profile with c= a. In the latter case kj = 0, and the foot of the per- 
pendicular from the virial focus Gy on the lift force is the Hamilton Centre C. 
We have in this case CP = a sin a, and the co-ordinates of C are (a sin? a, 
a sin a COs a), i.e, we have 
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THEOREM 6.—For the straight slit profile, the locus of the Hamilton Centre 
is a circle having as its diameter the line joining the foci of the meta- 
centric and virial parabolas. 


In the former case k, = (a? — c?)/2a and we easily get 


alll ko 
CP=asina+ = 


and the co-ordinates of C are given by 
Xo =k, +asin?a; ye= ky) cota+asinacosa (27) 


Elimination of a gives the locus of C as a cubic with symmetry about an 
axis perpendicular to the first axis, and touching the first axis at the centre 
of the profile. We can also trace this curve by using polar co-ordinates with 
PC =r, and 6 = 90° — a, getting 


r=acos 6 + ky sec 0 (28) 


which can be traced in a manner analogous to Rankine’s method of super- 
position of streamlines by using the two curves r, = a cos @ which is a circle 
and r, = ky sec 6 which is a straight line, and taking r = r, + rg. 


In the present case we can also find a minimum value for CP. Here 
CP = asina-+k,coseca and elementary analysis shows that this is a 
minimum for sin? a = k,/a and that the minimum is given by 


(PC) min = V 2 (a? — c?) (29) 
For the Karman-Trefftz profiles obtained by the transformation 


2+ _ (ts) 


yaa z’— Co 


(30) 


where ¢,, Co, are constants, and 7 is nearly equal to 2 (von Mises, 1941, 
p. 62), so that we can write n = 2 + «, where « is small, we have as in the 
case of the Joukowski profile y= 0, and 


ko = x (a? — c* cos 28). 


Here we can further write 


2 ——_ 
kya d {ar — 2 =U oo 2p) 


Hence ky > 0 if 


a e 
a a,” cos? B cos 2B < a’, since cos B = C2/ay 





B. S. MADHAVA Rao 


n?— 1 a? 
a,? <a’, or, —> 
3 a 


n? — | 


3 


Since a > a,, this condition can always be satisfied by suitably choosing G, 
in terms of which 4(m? — 1) = 1 + $e + 46%. Thus the same 
conclusions regarding the relative positions of C and P which were derived 
in the case of the Joukowski profile hold in this case also. 


Finally we come to the S-profiles given by von Mises (1941, p. 66). 
Confining ourselves to the case where there are only three singular points 
situated inside the circle at the points a, = c’,a.=c’e®, a; = — c"e¥, 
we derive the transformation 

, , 2 + oe"? e2t6 (c'c"e'@)?2 
zZa=2+ ra —_- — ia 
Zz a (31) 
which has the property that 8 = y, i.e., a profile having a centre of lift, or 
where the metacentric parabola reduces to a point. To find the sign of ky 
in this case, we notice that in the usual notation, B’G’? = c‘/c’*, where 
OG’ = c"/c’ and c’ angle (Fig. 5) x OG’ = twice angle x Oa,, since 


crett? — o'2 + c” 2210, 


Hence B’G’ = c?/c’. Since for the S-profile 8B = y, we shall have k, > 0 
if a >, i.e., a? > c*® or c’.B'G’ <a®. Now, we can choose aj, a2, a; and 
the circle of sufficient radius so as to pass through — c’, and include the 
singular points. Having made this choice, we can further alter it, preserving 
the original conditions, and so choose the centre M such that a?/c’> B’G’, i.e., 
such that k, > 0. Thus for an S-profile we can make the Hamilton Centre 
lie above the metacentre as in the case of the symmetric Joukowski profile. 
Hence for the S-profile the situation regarding stability is more favourable 
than in the Joukowski and Karman-Trefftz profiles. Although the situa- 
tion is the same in the cases of the symmetric and S-form profiles, we have 
this difference that in the former case 8 = y = 0, and the lift force is very 
small. 


In the above S-form profile we had 4) = 0 and k, > 0 but we can show 
that it is possible to find more general profiles of this type with h, > 0 and 
ky > 0, i.e., where the metacentric parabola opens upwards and C lies above 
P or the situation regarding stability is more favourable than even in the 
case of the S-form profile with a centre of lift. For this purpose let us write 
in (31) 6? = c"*e#4@ and let us take the centre of the circle, viz., M (Fig. 5) 
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as given by the complex number m = Xo + iyo such that a,a, passes through 
M, i.e., we take here 


mM = Xo + ivy = kb = ke"e’@ (32) 


where k is a real constant. The first axis is the line joining the singular point 
—c’ to the centre. Hence 





= , “= 2Vo (x + c’) 
enaits Xy +c” caeeciealts (X_ + c')® — ye* 


c”? sin 20 2x, ; 
c® +c’ cos 20 kc? rig =a using (32) 


tan 2y = 





tie , - > 
Hence the condition for hy = 0, viz., y = B, or tan y = tan B reduces after 


a slight simplification to the condition 


Xo" + Yo? + Xo (1 — kK) c’ = kc? (33) 


The equality gives the usual S-profile with the centre of lift. The first sign 


of inequality gives hg > 0. The last sign of inequality gives h)< 0 like the 
Joukowski profile. 


Further regarding ky 20, this is given by 


a? 2=c* cos 2(y — f) (34) 
This condition becomes rather complicated when expressed in a form ana- 
logous to (33). We therefore make a further simplifiction of (32) by taking 

Xy=Yo=me = (w > 0) (32’) 
With this condition (34) reduces after some simplification to 

ee 34’) 


and since » > 0, and both k and » are left arbitrary, a choice is possible such 
that the upper sign of the inequality in (34’) holds, i.e., ky > 0. Hence (33) 
and (34’) show that it is possible to choose a profile similar to the S-profile 
with hy > 0 and k, > 0. 





SUMMARY 


Problems of the virial associated with a simple wing profile are studied 
by considering the relative positions of the Hamilton Centre and metacentre. 
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In particular, questions of stability are emphasised. Several theorems are 
derived relating to the virial parabola, the Hamilton Centre and the axes 
of the profile. These theorems are specialised to the case of the Joukowski 
profiles, particularly the symmetric and straight slit cases. The S-form pro- 
files are generalised to the case where the metacentric parabola opens up- 
wards, and the Hamilton Centre is above the metacentre, which are the condi- 
tiens most favourable for stability. 
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ON THE INTENSITY AND POLARISATION OF THE 
LIGHT FROM THE SKY DURING TWILIGHT 


By J. V. DAVE AND K. R. RAMANATHAN, F.A.Sc. 


Received December 24, 1955 


1. INTRODUCTION 


THE object of this paper is to present the results of some new homogeneous 
measurements of the brightness of the clear zenith sky taken at Mt. Abu 
(24° 26’ N; 72° 43’ E) with a photon multiplier and six narrow band filters for 
different values of solar depression from 0° to 20° together with some 
measurements of polarisation, discuss the results in the light of our present 
knowledge of the upper atmosphere and attempt on broad lines an 
explanation of the optical phenomena observed during the transition 
between day and night. 


2. EXPERIMENTAL ARRANGEMENTS 


An R.C.A. photomultiplier 931-A, selected for small dark current, was 
used with an A.C. operated D.C. power-supply stabilised by a degenerative 
type of electronic regulation. The photomultiplier output current was ampli- 
fied by a balanced-bridge stabilised D.C. amplifier. The sensitivity of the 
amplifier could be changed successively by factors of about 10 by introducing 
different resistances varying from 10* to 10° ohms in its grid circuit. The 
exact scaling factor from one stage to the next was determined by illuminat- 
ing the photocathode by a beam of constant intensity and noting the reading 
in each successive stage. No attempt was made to measure the absolute 
values of the intensities; only relative intensities were determined. For 
values of output current less than 400 microamperes, the current was found 
to be proportional to the intensity of light falling on the photocathode. For 
higher outputs, the correction due to the non-linearity of the variation of 
output against input voltage was measured and applied. The effect of vary- 
ing the orientation of the plane of polarisation of the incident beam on the 
sensitivity of the photomultiplier was tested and found to be negligible. 


Suitable optical filters were used to transmit different parts of the spec- 
- trum. The components and characteristics of the filters are shown in Table I. 
The figures take into account the spectral sensitivity of the photomultiplier. 
It will be noted that with the green, yellow and red filters, the effective trans- 
missions are about 200A centred approximately at the respective lines of 
the airglow spectrum. 
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TABLE I 


Components and effective transmissions of optical filters when used 
with RCA Photomultiplier 931 A 





Brief name 


Components 
for filter 


Half band-width 
of transmission 


Maximum 
transmission 





Chance OX | 
Chance OV | 
Chance OY 18 and OB 10 


3300-3850 
3640-4050 
4350-4800 


3700-32% 
3850-44 
4600-94 


LF. plus ON 16 plus Plastic Yellow 5530-5750 5600-34% 
5870-2% 
6300-19 


4600-40% 


1.F. plus Chance OY | 5800-6000 


LF. plus Plastic Red 6150-6400 


Polaroid 4000-5500 





I. F. stands for interference filter obtained from Messrs. Barr and Stroud. 


For the measurement of polarisation, a polaroid of type J was placed 


in the path of the radiation so as to measure (1) the intensity I, of the compo- 
nent perpendicular to the sun’s meridian and (2) the intensity I, of the 
component parallel to it. Both the intensities were measured within a time 
interval of 5 seconds. The percentage polarisation was calculated from the 
formula, 


3. RESULTS 


(a) Intensity —Observations were made only on cloudless clear skies. 
Only one narrow-band filter was used on any particular morning or evening. 
Meter readings were noted at regular intervals of half a minute except during 
late twilight when the rate of change of intensity was very small. Frequent 
check was made on the zero reading and photomultiplier voltage. After 
applying linearity and zero shift corrections, the corrected meter readings 
were multiplied by appropriate scaling factors so that the final values were 
proportional to the intensity of the light falling on the photocathode. 


The readings taken on different clear days with the same filter were 
found to agree well when the photomultiplier voltage was kept very steady. 
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The actual readings for the same solar depression were found to differ in 
certain ranges on some days. This might have been due to srhall changes 
in photomultiplier voltage, changes in the scattering and absorption of light 
by dust and haze in the atmosphere and changes in the emission of the air- 
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Fic. 1. Intensity of the light from the zenith sky in evening twilight in defined regions 
of the spectrum. 
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glow lines. For comparing the course of the variation on different days and 
in different wave-lengths, the readings were reduced to a common standard 
by multiplying the values on a particular morning or evening by a factor 
which would make the readings on all the occasions have the same value at 
a solar depression of 5°. The means of the intensities observed on different 
days with the same filter and for the same solar depressions, were calculated. 


The average values of the observed intensities of the zenith sky in arbi- 
trary units during evening twilight in different portions of the spectrum are 
plotted against solar depression in Fig. 1., The intensity scale is logarithmic. 
The following are the main features of zenith sky brightness during twilight: 


(1) The curves for the different wave-lengths from the near ultra- 
violet to red are generally similar from D = 0 to D = 12°, but in the shorter 
wave-lengths, the rate of decrease of intensity with increasing solar depression 
becomes markedly less when D is about 8°. 
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Fic. 2._ Comparison of twilight ‘zenith sky intensities as measured i in emacnase places. 
a ar ial ies nis + ae mee 
~ GGALy +4400 A. Ashburn (Cactus Peak, California). 
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Fig. 3. Intensity of parallel and perpendicular components of zenith sky light during 
morning twilight (4000-5700 A). 


A. Parallel component. B. Perpendicular component. 
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(2) While the intensities of the ultraviolet, violet and blue flatten out 
when the sun’s depression is 16° ta 18°, the intensities of the green, yellow and 


red which include bright lines in the night airglow flatten out when D is 13° 
to 14° at a higher level of relative intensity. 


(3) In the green and the red, there were observed on many days a reduc- 
tion in the rate of fall of intensity when D was about 12° suggesting a tempo- 


rary brightening near the boundary of earth’s shadow. This requires fur- 
ther examination. 


For comparison, the values of zenith sky intensities observed during 
twilight by other observers are plotted side by side with the Abu values, 
Ashburn’s measurements! in California were made with a phototube and 
filters and come closest to Abu values. The observations of Chiplonkar and 
Ranade? were made at Sinhagad near Poona with a visual photometer includ- 
ing the oxygen green line. Ljunghall’s* curve is composite being made up 
of his own observations at Helwan with a photocell and filter and of the old 
visual observations of Dufay and others. 


(b) Polarisation —Observations of I, and I, for studying the polarisa- 
tion of the sky were taken on a number of mornings and evenings: (1) with 
the polaroid alone (which transmits the violet, the blue and the green includ- 


ing an appreciable fraction of the green line 5577 A) and (2) with polaroid 
and blue filter. 


In Fig. 3 are plotted the values of I, and I, for different solar depressions 
during morning twilight with polaroid alone. The corresponding curves 
for evening twilight taken with polaroid and blue filter for the spectral region 
4350-4850 A, are essentially similar. 


The values of the percentage polarisation as calculated from the values 
of I,, and I, are plotted in Fig. 4. Curves B and C represent the values 
obtained in the effective spectral region 4350-4850 A for evening and morn- 


ing twilight respectively, while the curves D and E are for the spectral region 
4000-5700 A. 


The curves showing zenith sky polarisation against solar depression are 
in general agreement with those obtained by Robley* visually at Pic du Midi. 
His observations extended to 16° and are also -shown in Fig. 4, curve A. 
From Fig. 4, the following points are clear. When observations were taken 
with a blue filter, the zenith sky polarisation remained more or less steady 
at about 70% till D was 5°, and then showed a pronounced fall till D was 
84°. Between 84° and 12° solar depression, the polarisation again showed a 
plateau at 45% to 50%, after which it declined rapidly so that when D was 
greater than 18°, the zenith sky was practically unpolarised. 
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Fic. 4, Percentage polarisation of the light from the zenith sky during twilight. 
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Robley’s curve shows similar variations, but the initial polarisation at 


low values of solar depression is higher, due no doubt to the clearer skies 
over Pic du Midi. 


The curves obtained with the polaroid alone without the blue filter show 
similar features but the changes of polarisation are less sharp. 


The polarisations at points 30° from the zenith in the vertical plane 
through the sun were measured both towards the sun and away from the sun 
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Fic. 5. Percentage polarisation of skylight in morning twilight (4350-4850), 
A. 30° from the zenith away from the sun. 
B. At zenith. 
C, 30° from the zenith towards the sun, 








Intensity and Polarisation of Light from the Sky during: Twilight 75 


on different days during morning twilight. The values of the polarisation 
calculated from the observed values of I;, and Ip are plotted in Fig. 5. 
Curves B, A and C represent the variations of polarisation with solar depres- 
sion at zenith, 30° from the zenith towards the sun and at 30° from the zenith 
on the anti-sunside respectively. 


It will be observed that the changes in sky polarisation 30° away from 
the zenith but towards the sun are sharper than those at the same angular 
distance from the zenith on the anti-sunside. The corresponding changes 
occur as may be expected, first on the anti-sunside, then at the zenith, and 
lastly on the sunside; the flattening out of the polarisation curve occurs 
at 9-5° on the sunside, at 9° at the zenith and at 8-5° on the anti-sunside. 


To explain the humps in the polarisation curves requires further ana- 
lysis of secondary scattered light. 


4. DISCUSSION 


When the sun is a few degrees below the horizon, the upper part of the 
atmosphere is directly illuminated by the sun’s rays, while the lower part 
lies in the shadow of the earth. The shadow region is however illuminated 
by scattered radiation from the sky all round, particularly from the hemi- 
sphere towards the sun. The light from the zenith sky is therefore made up 


of the primary scattered light from the upper part of the atmosphere and 
secondary (and multiple) scattered light. Owing to the greater air density 
in the lower part of the atmosphere, however, the main secondary scattered 
light would come from layers lower down in the atmosphere than the primary 
scattered light. As a consequence, as the sun goes down, the primary 


scattered light would decrease at a much faster rate than the secondary 
scattered light. 


The observed variation of the brightness of the zenith sky with the de- 
pression D of the sun below the horizon agrees with the view that the sky 
illumination during twilight is only in part due to primary scattered light P. 
If it were due mainly to P, the variation of the sky brightness with D would 
be similar to the variation of the amount of air above the lower boundary 
of the directly illuminated upper atmosphere, that is, of the pressure at the 
shadow boundary. In Fig. 6 are drawn the curves of brightness RR and BB, 
of the zenith sky observed at Abu in the red (6300 A) and the blue (4600 A) 
and also the curve of pressure PP at the shadow limit. The height of the 
earth’s shadow above the ground (for 4500 A) corresponding to different 
values of D are taken from Ljunghall’s paper. Ljunghall assumes that above 
6km., the attenuation of light in the atmosphere is due only to molecular 
scattering while at lower levels, the coefficient increases, becoming at ground 
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level 2-2-5 times the value for molecular scattering. The refraction due to $0 
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Fic. 6. Comparison of the observed zenith sky intensity with pressure at the base 

of the illuminated layer. 

RR Observed intensity at 6300A. 

BB Observed intensity at 4600 A. 

PP Pressure at the base of illuminated layer. 

Height of base as given by Ljunghall. 

SS_ Integrated illumination at ground level due to scattered light from the 

whole sky during twilight, i 
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solar depression of about 6°, they deviate more and more from it at greater 
values of D. 


At solar depressions of 12° and 15°, the values of BB are more than 
a thousand times greater. Many years ago, it was suggested by Hulburt® 
that secondary scattering plays an important part in twilight illumination. 
He showed that when the solar depression exceeds 8°, the light scattered from 
the zenith sky from the lower part of the atmosphere lying within the shadow 
of the earth but illuminated secondarily by the twilight glow is comparable 
in amount to the light scattered by the upper part of the atmosphere directly 
lit by the sun’s rays. We now find that the discrepancy is enormous; can 
all this be explained by multiple scattering ? From observations now avail- 
able, it is possible to make an estimate of the relative contributions of prim- 
ary and multiply scattered light to the brightness of the zenith sky during 
twilight. 

Koomen® and others have given values of sky brightness measured in 
different azimuths and altitudes at Sacramento Peak in New Mexico (2800 m. 
above sea-level) for different values of sun’s depression down to 15° below 
horizon. From these measurements, one can by numerical integration 
estimate the illumination due to scattered light from the whole sky at the 
point of observation on the earth’s surface. This has been done and the 
following are the results. 








Sun’s depression Illumination 
below due to skylight 
horizon at ground 
0° 180 candles/ft.? 
Ki 22°6 
6° 0-76 
‘ 9° 0-031 
12° 0-0017 
15° 0-00027 





The values in column 2 are plotted in Fig. 6 as curve SS. 


We may assume as a first approximation that the scattered light incident 
after sunset or before sunrise on a vertical column of the atmosphere above 
the observer will vary similarly to the integrated illumination at ground as 
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given in column 2 above. PP and SS will then represent respectively the 
relative variations of primary and secondary scattered light from the atmo- 
sphere during twilight. It is known that when the sun is on the horizon, 
i.e., at D = 0, the secondary scattered light from the zenith for 4500 A is 
nearly 0-25 of the primary scattered light, but we do not know how this 
ratio will change when D increases. In what proportions the primary and 
secondary should be combined at different times after sunset is a matter for 
investigation. 
5. POLARISATION OF LIGHT FROM ZENITH SKY 
DURING TWILIGHT 


The view that secondary scattering is the main factor responsible for 
extending the twilight period is borne out by the measurements of polarisa- 
tion of the sky during twilight. Robley first showed from his observations 
at Pic du Midi that the polarisation of the clear zenith sky during twilight 
fell rapidly when D increased to more than 6° but that it later attained 
a steady value of 50% to 55% when the sun’s depression was between 
9° and 12°. These observations have been corroborated by our Abu 
observations (Fig. 4). When the secondary scattering also becomes in- 
significant, the light from the zenith sky is only scattered starlight and light 
from the airglow and is practically unpolarised. This happens when D is 
16° to 18° for blue light. Robley found that with red light for which the 
secondary scattering would be much smaller, the first dip in polarisation 
and the plateau between D = 9° and D = 12° were inconspicuous. It may 
be observed that the sky illumination during evening twilight, especially in 
its later half, is markedly anisotropic, the light coming from the direction of 
the sun being much greater than from other directions. This light, when 
secondarily scattered downward by the atmosphere above the observer, will 
naturally have a large polarisation. 


Quantitative checking of these ideas with observation requires the 
calculation of the intensity and polarisation of secondary scattered light 
by the spherical atmosphere when the sun is below the horizon. The calcu- 
lations have been made and will be the subject of communication of Part II. 
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Ir has been reported earlier that 7-methyl ethers of polyhydroxy isoflavones 
could be. prepared both by partial demethylation'~* as well as partial methyla- 
tion.4 A method which has been successfully employed in the preparation 
of 5-methyl ethers depends on the partial benzylation of the other reactive 
hydroxyls, methylation of the 5-hydroxy group and final debenzylation.? 
Recently benzoylation® has also been employed in connection with a synthesis 
of muningin. It has, however, been observed in this laboratory that a more 
convenient method is to effect partial acetylation of the reactive hydroxyls, 
methylation in the 5-position and final deacetylation.® 


Among naturally occurring flavonoids, muningin (XJ) is unique in hav- 
ing a methoxy group in the 5-position and having at te same time more 
reactive hydroxy groups free in other positions. Further, it is a 5: 7- 
dimethoxy compound and hence a combination of methods was needed for 
its synthesis and for the synthesis of analogous compounds. As the first 
step the preparation of the 7-methyl ether (III) of 5:6: 7-trihydroxy iso- 
flavone (I) was attempted. It was observed that partial demethylation of 
5: 6: 7-trimethoxy isoflavone with hydriodic acid did not proceed as smoothly 
as with members of the 5: 7-dimethoxy type and the yield of the required 
product was low. As an alternative, partial methylation of 5: 6: 7-tri- 
hydroxy isoflavone (I) was attempted using the acetone, potassium carbonate 
method and 1 mole of dimethyl sulphate. The major product was 6: 7- 
dimethoxy-5-hydroxy isoflavone (II) and the monomethyl ether (III) could 
not be successfully obtained. Even in the simpler case of 2-methyl-6: 7- 
dihydroxy isoflavone (IV) partial methylation could not be effected by this 
method and only dimethyl ether (V) could be obtained. Obviously the 
difference in the reactivity of the 6 and 7-hydroxyl groups is not appreciable 
enough for this method. 


On the other hand by avoiding potassium carbonate and using sodium 
bicarbonate instead, and carrying out the methylation over a longer period, 
79 
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the 7-monomethyl ether (VI) could be obtained without serious difficulty. 
The bicarbonate method was used earlier’ for a similar case in the coumarin 
series, that is, partial benzylation of 6: 7-dihydroxy coumarin. Extending 
the method, 7-methoxy-5: 6-dihydroxy isoflavone (III) could be conveniently 
prepared, and when it was subjected to partial acetylation of the 6-hydroxy 
group, methylation of 5-hydroxyl and final deacetylation it yielded 5: 7- 
dimethoxy-6-hydroxy isoflavone (VII) whose properties were the same as those 
described earlier® except the ferric reaction. Following the same procedure 
and starting from 5:6: 7:4’-tetrahydroxy isoflavone (VIII) the 7-methyl 
ether (IX), and its diacetate (X) were obtained and final methylation and 
hydrolysis gave — (XT) (see also ref. 5 and ing 
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5:6: 7-trimethoxy and 5:6:7:4’-tetramethoxy isoflavones were pre- 
pared from the appropriate phenyl benzyl ketones by ring closure using ethyl 
formate. Yields upto 20% were reported earlier in the isoflavone con- 
densation of 2-hydroxy-4: 5:6: 4’-tetramethoxy phenyl benzyl ketone by 
Krishnamurthy and Seshadri® and King et al.’°; it has now been found that 
if the reaction mixture is kept for 96 hours instead of the usual 24 hours the 
yields are as good as 55%. 5:6: 7-Trihydroxy isoflavone and 5: 6:7: 4’- 


tetrahydroxy isoflavone were obtained in very good yield by demethylation 
with hydriodic acid, 
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1-Methoxy-5 : 6-dihydroxy isoflavone (III) 


5:6: 7-Trihydroxy isoflavone (2 g.) was dissolved in acetone (100 c.c.), 
methyl sulphate (0-7 c.c.; 1 mole) and sodium bicarbonate (5 g.) were added 
and the mixture refluxed for 24 hours. Acetone was distilled off and the 
residue treated carefully with dilute hydrochloric acid till all the bicarbonate 
was decomposed. After filtering, washing and drying, the product crystallised 
from alcohol as fine silky needles melting at 222-24°. A mixed melting 
point with an authentic sample prepared by Aghoramurthy, Venkatasubra- 
manian and Seshadri!’ was undepressed. Yield 0-7 g. 


7-Methoxy-6-acetoxy-5S-hydroxy isoflavone 


The above methyl ether (1 g.) was dissolved in pyridine (8 c.c.) and acetic 
anhydride (0-5 c.c.; 1-5 moles) added. The solution was stirred well for 7 
minutes. During this process a solid started separating out. Pieces of 
ice and water were added and the product was filtered, washed well and 
dried. It crystallised from ethyl acetate as colourless rectangular plates 
melting at 184-85°. Yield 0-8g. It gave a red colour with alcoholic ferric 
chloride (Found: C, 66-5; H, 4-6; C,sH,,O6¢ requires C, 66-3; H, 4-3%). 


5: 7-Dimethoxy-6-acetoxy isoflavone 


The above partial acetate (0-5 g.) was methylated in acetone solution 
using methyl sulphate (0-4 c.c., excess) and potassium carbonate (2 g.) for 
40 hours. The product crystallised from alcohol as colourless thick rectangu- 
lar tablets melting at 205-07°. Yield 0-3g. (Found: C, 67-7; H, 5-1; 
Ci9H,6O¢ requires C, 67-1; H, 4-7%). 


5: 7-Dimethoxy-6-hydroxy isoflavone (VII) 


The above acetate of the dimethyl ether (0-2 g.) was placed in a small 
test-tube immersed in freezing mixture. Ice cold sulphuric acid (2 c.c.) was 
added dropwise during the course of about five minutes, and the tube was 
then transferred to another beaker containing ice and water. It was stirred 
and kept at this temperature for 10 minutes. Addition of ice caused the 
separation of the deacetylated product as a colourless solid which crystallised 
from dilute alcohol as colourless thick rectangular tablets and prisms 
melting at 184-86°. It was freely soluble in dilute alkali but gave no colour 
with ferric chloride (Found: C, 68-1; H, 5-2; C,,H,,0, requires C, 68-4; 
H, 4-7%). 


7-Methoxy-5: 6: 4'-trihydroxy isoflavone (IX) 


5: 6:7: 4'-Tetrahydroxy isoflavone (2g.) in acetone solution was 
methylated with methyl sulphate (0-7 c.c.; 1 mole) and sodium bicarbonate 
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(5 g.). The product crystallised from acetic acid as small prisms melting at 


263-65° (decomposition) with shrinking at 258°. Venkataraman et al.’ 


reported the same melting point. Yield 0-6g. (Found: C, 64-0; H, 4-3; 
C,6H,. Og requires C, 64-0; H, 4-0%). 


7-Methoxy-6 : 4'-diacetoxy-5-hydroxy isoflavone (X) 


The above mono-methyl ether (0-5 g.) was partially acetylated with 
acetic anhydride (0-4.c.c.; 2-5 moles) in pyridine (5 c.c.) solution as before. 
The product crystallised from ethyl acetate as colourless stout rectangular 
prisms melting at 206-08°. It gave a red colour with ferric chloride; yield 
0-4g. (Found: C, 62-5; H, 4-5; C.pH,sOs requires C, 62-5; H, 4-2%). 


5: 7-Dimethoxy-6: 4'-dihydroxy isoflavone (muningin) (XT) 


The above compound (0-5 g.) in acetone solution was refluxed with 
dimethyl sulphate (0-5c.c.; excess) and potassium carbonate (4g.) for 
40 hours. The product separated from alcohol as colourless needles melting at 


230-32° agreeing with the melting point recorded for muningin acetate by 
King et al.° 


The above acetate (0-2g.) was deacetylated using cold concentrated 
sulphuric acid. The product gave all the reactions for muningin as de- 
scribed by King et al."° It crystallised from dioxane as colourless rectangular 
tablets melting with decomposition at 285°. Mixed melting point with the 
natural sample was undepressed (Found: C, 64-5; H, 4-8; C,;H,4O¢ requires 
C, 65-0; H, 4-5%). . 


SUMMARY 


Of the two methods of obtaining 7-methyl ethers of 5:6: 7-hydroxy 
isoflavones partial methylation employing sodium bicarbonate is more 
convenient. Thus 7-methoxy-5: 6-dihydroxy and 7-methoxy-5: 6: 4’-tri- 
hydroxy isoflavones have been prepared. These could be partially acetylated, 
methylated in the 5-position and finally deacetylated to yield 5: 7-dimethoxy- 
6-hydroxy isoflavone and muningin respectively. 
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ABSTRACT 
A symmetric oscillatory Poisson process is defined and its stochastic 
features studied. The process represented by the symbolic integral of 
this oscillatory Poisson process is then discussed in detail. The results 
obtained are applied to the well-known stochastic problem of multiple 
scattering of charged particles in their passage through matter. 


INTRODUCTION 


In this paper we shall study processes which are associated with what we 
shall call a symmetric oscillatory Poisson process (hereinafter referred to as 
S.O.P.) defined as follows. Events occur in a Poisson manner, along the 
one-dimensional ‘ t’ axis, i.e., if n(t) is the stochastic variable representing 
the number of events in the interval (0, ¢), then the probability that n(t) =n 
is given by = (n, t) = e~** (At)"/n!, where A is the parameter of the Poisson 
distribution and Adt, the probability that an event occurs between ¢ and 
t+ dt. If we associate with the i-th event a stochastic variable x; which 
assumes the values + 1 and — | with probability p and (1 — p) respectively 


i=n (0) a js e 
then m(t) = & x; is a stochastic variable which can assume integer values 
i=0 


from — co to + co and represents an oscillatory Poisson process (O.P.). 
If p = 4 we call it, for obvious reasons, an S.O.P. The typical trajectory, 
that is, the realised curve of n (r) in the interval o to t, r representing a typical 
point on the ¢-axis, consists of lines parallel to the ¢-axis with jumps of unit 
magnitude (positive and negative) at points where the Poisson events occur. 


We shall here deal with the following: 
(1) The probability distribution of m (t), 
(2) Recurrence and first passage times of m, 


(3) The processes represented by the symbolic integrals 
Fy m(r) dr and f m(t— dr, 
0 0 


(4) Application of these results to the problem of multiple scattering, 
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(5) Generalisation of m to the case m = J a; where a; is a stochastic 


4 
variable associated with the i-th event, all the a;’s having the 
same probability frequency function. This problem is shown 
to be equivalent to that of Snyder and Scott (1949). 


(6) Multiple scattering in multiplicative processes. 


The object of the paper is twofold: (1) First is to explain, in physical 
terms, the importance of the concept of trajectory in the computation of 
first passage and recurrence times. The authors have not seen a systematic 
treatment of this concept and its relevance to the problems of first passage 
and recurrence from a physical point of view and have assumed that a summary 
of the present state of knowledge is contained in Bartlett’s recent book (1955) 
on stochastic processes. The sections on recurrence times in this paper 
amount to an extension of the ideas in Bartlett’s book but with additional 
emphasis on the concept of the trajectory of a stochastic process and conse- 
quently on the necessity of a distinction that has to be made between two 
“equivalent ’’ processes (to be explained presently) when considering first 
passage and recurrence times. This aspect has not been stressed in Bartlett’s 
book. The phenomenological difficulties relating to first passage and 
recurrence problems have been classified and discussed in a recent contri- 
bution by one of us (R, 1955) and the general considerations of that paper 
are applied to the solution of the problems here. 


(2) The second object is to consider the problem of multiple scattering, 
as an integral of a basic random process and suggest an approximation which 
helps us to obtain exact solutions. The last section on multiple scattering 
and multiplicative processes deals with the beginnings of one of the most 
vexed problems of cascade theory in cosmic radiation, i.e., the lateral spread 
of showers. 

1. THE PROBABILITY DISTRIBUTION OF m(t) FOR AN O.P. 

The standard procedure is to express 7 (m, t + dt)* the probability that 
m(t + dt) =m in terms of 7(m, t) that probability that m(t) = m. 

The forward differential equation of the process is then given by 

om (m, t) _ 


rT; — Ax(m, t) + A{pr(m — 1, 1) + (1 — p)z(m + 1, D}, 
a(m,0)=0, (m+0) 
a (0, 0) = 1. (1) 





* Throughout this paper, x denotes any probability frequency function; distinction 
between two functions will be apparent from the context. Where Fourier transforms are 
defined for these functions, the same convention is adopted. 


A3 
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It is not easy to solve the above equation by simple iteration, since 
m(m, t) is expressed in terms of both 7(m-+ 1, t) and 7(m—1, 12). The 
validity of defining a generating function 


G (u, t) = L7(m, t)u™ (2) 


depends upon the uniform convergence of the series with positive and nega- 
tive powers of wu in a domain including u=1. To avoid these difficulties 
we adopt the following argument. 


The above process can be considered as a superposition of two Poisson 
processes, having the parameters Ap and A(l — p) but while one of these 


has a value + 1 associated with each event the other has the value — 1 
assigned to each event. 


Using this argument we immediately obtain the frequency function 
a (m, t) as 


r=Co 


m(m, t) = BD ee er e- (1-p)t oe (3) 


r=0 Or —m 


In the summation, the lower limit of r is zero if m is positive and — m if m 
is negative. In particular if we assume p = 4, i.e., m(t) is an S.O.P., 


a (m, t) = e-™ Im (At), (4) 
where I, is the modified Bessel’s function of the first kind of order m. Note 


that a (m, t) = 7(— m, f), i.e., the distribution is symmetric. about the value 


m= 0. Then « {m(t)} = 0* and « {m?(t)} = At, the same as the mean square 
deviation of n(t) with parameter A. 


2. RECURRENCE AND First PASSAGE TIMES OF m 


We now ask for the probability F (m, t) dt that the variable m (7) reaches 
the value m for the first time between ¢ and ¢ + dt. For any m <0, F (m, t) 
can be obtained by writing 7 (m, ¢) in terms of F (m, t) which is easy for any 
Markoff process. By simple arguments, we have 


a(m, t) = fF(m, 1) 7 (m| m, t — 7) dr (5) 


where 7(m|\m, t) is the probability of having a state m at t, given that the 
state at 0 is m. In our case as the variable m(r) is purely additive, this pro- 
bability is identical with 7 (0, tf) and hence we can immediately obtain a solu- 
tion for F (m, t) by resorting to the Laplace transform technique. Defining 


+ Throughout this paper e denotes the expectation valuc, 
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p(m, s) and f(m, s), as the Laplace transforms of 7 (m, t) and F (m, t) res- 








pectively, we have : 

p (m, s) = f (m, 5) p (0, 8) (6) 
or 

f(m, s) = p(m, s)/p (0, s) = f(— ™, 5). (7) 
Now from (1) 

_ &+- yop” 

P(m, 8) = a (8) 
Hence 

f(m, s) =f(— m, ) = SEAVER IE ) 
which yields 


F (m, t) = F(— m, t) = cd (10) 
This expression is valid for all m except m= 0. In the latter case a 
first passage is said to occur only if the state has changed from o at some 
point between o and ¢ and then returned to o for the first time at ¢ and this 
is more appropriately called a “‘ recurrence”. Continuance in the state o 
from o to t + dt does not constitute a recurrence between ¢ and t+ dt. On 
the basis of this definition, we find that the state o may be realised at ¢ either 
by a first passage to o at some t < ¢ and a subsequent realisation of the state 
o at t or directly by a continuance of the state o throughout the interval o to t. 
We are thus led to the equation 


(0,1) = fF (0,1) (0,t—2)dr-+e™ (1) 


so that 





P (0,8) =f(0, 8) P(0,8) + ya 
On using (5) we obtain 


AW 6/ie tie — 9 
fo,s) = s+ vist +)? — (12) 





and hence 


t 
F(o,t) = de f LO dr. 
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F (0, t) dt gives the probability that the state o recurs between ¢ and 
t + dt and F (0, t) therefore gives the probability frequency function of the 
** Recurrence Time ”’. 


If we try to obtain the mean value of the first passage time by multiplying 
the R.H.S. of (10) by ¢ and integrating with respect to ¢ or directly from 
f(m, s), we are led to the interesting result that the first, and consequently 
higher, moments of the first passage times are infinite for any m. 


t 
3. THE PROCESS REPRESENTED BY THE SYMBOLIC INTEGRAL / mz) dr 
0 


t 
To study the process symbolically represented by the integral / m (rt) dr 


we adopt the phenomenological interpretation of integrals of random func- 
tions given recently by one of us (R., 19555, c). A typical realised tra- 
jectory of m (7) is given in Fig. (1). 












































Te 








wet 
Fic. 1. Typical trajectory of m (z) in the interval (0, 1). 

For a given realisation of random points along the ¢ axis in the interval 

o to t, the corresponding realised value of y(7), which we shall denote by 
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y®(z), is equal to { m®(z') dr’ where m® (r’) is the realised value of m (7'). 


y®(z) represents the area enclosed by the realised curve of the stochastic 
variable m(7’) as 7’ varies from o to 7. In computing the area, the area 
below the 7 axis has to be taken as negative. Thus we plot the realised curve 
of y (7) in the interval o to t; see Fig. 2. 














T 
Fic. 2. Typical trajectory of y (1) = J m(z’) dz’ in the interval (0, 1). 


We note that it is continuous and its slope at any point 7 is the realised 
value of m® (7) at that point. y therefore is a continuous stochastic variate 
and we can define the probability 7(y, t) dy that y(t) lies between y and 
y+dy. The obvious method of obtaining the equation for 7(y, ft) is to 
study the variation of the area referred to above by increasing ¢ to t + df. 
Since the increase in area is m(t) dt we note that the process y(t) is non- 
Markovian and hence we are constrained to introduce the function 7 (m, y) dy, 
the joint probability that m(t) =m and y(t) lies between y and y + dy. 
The process defined by the pair m(t) and y(t) is Markovian and by simple 
arguments we obtain ; 
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om (ms) _ _ v(m, ys 1) + 5 fr (m — 1, y, 1) + (m+ 1, y, D} 
dm (m, y, t) 
— sea * era ’ (14) 
and 
+co 
7 (y, t) = En(m,y, 1). (15) 


It is difficult to solve this equation in view of the occurrence of m and 
to avoid this complication we take recourse to the method of regeneration 
points which yields an integral equation for 7(y, t) directly. 


r= { F(t 0-C-9,t- +70 4+0—9t- 9} 


dr+8(ye™. (16) 
This equation is obtained by using the now familiar arguments of the regene- 
ration point method that the first Poisson point occurs between 7 and r+ + dr 
with probability e~*7 Adr and the realised value of y when events are counted 
from 7 differs from that of y when counted from t=O by +(t—7) 
according as we assign + 1 or — | to the Poisson point at r. Differentiating 


the above equation we obtain what is known as the backward differential 
equation of the p.f.f. of y(t): 


e 
U9) dQ, )+5(7O-b)+e0+40. 


This equation is easier to solve. The equivalence of (14) and (17) is difficult 
to prove. But therein lies the advantage of dealing with equation (17)! 


Defining the Fourier transform of 7(y, t) as 
P(u, t) = [7 (y, t) ei dy, (18) 


equation (16) reduces to 


t 


Plu th= fer Zfetuen + etut- put — Ndr tem. (9 


0 


Differentiating this with respect to t we obtain 


PC 1) — — rp (ut) + 3 fet + p(w, 1), (20) 















A Symmetric Oscillatory Poisson Process 




















Hence 


p(u, t) = exp i- Mt + = (e ut — omy : (21) 











) Inversion of this expression yields 
2 x n — 2r)t}"-1 
a(y,t) =e et 5(y) + ert Tt 2 (* 1c prise = 
}) 
FI H(y +n—?2rt). (22) 
. The mean value of y of course is zero. Its variance may be obtained 
directly from (21); we find o? = Ar?/3. 
t 
3 (a). THE PROCESS REPRESENTED BY THE SYMBOLIC INTEGRAL f m (t—7) .7 
On examining equation (17), we observe that it can be interpreted a» the 
6) forward differential equation of another stochastic process which unlike 
; f a (7) dr is a quasi-Markovian (the significance of the qualification ‘‘ quasi ”’ 
t ry 
od will be explained presently) basic random process characterised by discrete 
7) jumps. For if we define a stochastic process by attributing a random 
ng variable x; to a point occurring at 7; where x; can take the values + 7; with 
al equal probability 4, then y* (t) = 2 x; defines a stochastic process progressing 
with t whose forward differential equation is given by (17). Its typical 
trajectory is given in Fig. 3. 
7) According to the theory of equivalent processes developed by one of us 
tt (R.), y*(t) can be represented by the symbolic integral 
l 
t 
1)! y*t() = fm (t — r) dr (23) 
0 
§) i.e., the integral is associated with the inverse trajectory of m(7). The con- 


cept of an inverse trajectory and the equivalent processes arising therefrom 
has been fully discussed by one of us in a series of three papers (R., 1955 5, c, d) 
and no purpose would be served here by a repetition of that discussion. 
Suffice it to say, that y*(t) and y(t) are two distinct processes, which pro- 
9) gress with ¢ in strikingly different ways as is revealed by the nature of their 
trajectories. But their probability distributions are identical at every t. 
y(t) is non-Markovian as has been explained in the derivation of (14) and 
(15). y*(t) is quasi-Markovian (as distinguished from Markovian), i.e., 
20) the state of y* (t) at t + dt can be predicted, if we know the state at t provided 
we assume we also know the origin of the process lies at t = 0, an information 
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y(t 




















%0 Tt 
T 
Fic. 3. Typical trajectory of y* (1) = f m(t — 1’) dz’ in the interval (0, #). The ordi- 
1 
nates have been scaled down for convenience of representation. 


unnecessary in Markovian processes. Otherwise, y* (t) has all the character- 
istics of the Markovian process in that no information between 0 and f¢ is 
needed, provided the state at t is known. It is obvious from the definition 
of y*(t) = XY x that its backward differential equation is just the forward 


differential equation of y(t). Hence the solution for 7 (y, ft) gives the distri- 
bution function of y*(t) also. 


First Passage and Recurrence Times of y(t) and y* (t) 


The above distinction between y(t) and y*(t) may seem academic if 
we are interested only in the distribution function. But as has been stressed 
by one of us, the distinction is fundamental, since the integral of a y (t) process 
is different from that of a y*(t) process and they are not even equivalent. 
The distinction is as vital when we consider first passage and recurrence 
times of the two processes. We shall consider the problems separately. 


y (t).—Since the trajectory of y(t) is continuou$ as shown by R., we 
have to define a function F (y, t), where F (y, t) dt represents the probability 
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that y(t) crosses the value y, between ¢ and ¢ + dt, or in other words, will 
assume exactly the value y, sometime between ¢ and t + dt. / F (y, t) dt = 1 
0 


for all values of y. The problem of determining F (y, ¢) is fraught with very 
great difficulties which have not been adequately resolved till now. But R. 
has given a simple expression for the mean number of recurrences in an 
interval o to t. We say that y recurs between ¢ and ¢ + dt when y(t) crosses 
y in the interval dt if it has crossed the same value earlier. It is of course 
assumed that continuance in the value y (this happens only when m is o in 
some interval) does not constitute crossing. As regards the problem of 
recurrence, we can only obtain, according to R., the mean number of re- 
currences in the interval (0, f). 


If the stochastic variable representing the number of recurrences of the 
state y is denoted by R (y, ¢), its mean value is given by 


co 


{R(y, D} = Emm (mm, y, 1). (24) 


Denoting the interval between the n-th and (” + 1)-th recurrences by Ty, we 
note that T,, is a stochastic variate and the distribution functions of the 
various T,’s are different for different values of n. So there is no meaning in 
speaking in general of a mean recurrence time. When the mean number of 
recurrences is very large, t/<{R (y; 1)} can be taken to be roughly an estimate 
of the mean recurrence time. 


y* (t).—Though trajectory of y*(t) is characterised by discrete jumps, 
as the magnitude of the jumps is a continuous variate, y*(t) is also a conti- 
nuous variate. It represents, of course, a basic random process. R. has 
shown that in the case of such processes we have to define F (y, t) dydt, as the 
probability that y (t) jumps into a value between y and y + dy, in the interval 
dt for the first time. Note the necessity. of attaching an infinitesimal dy to 
the function F, to give it the significance of a probability magnitude. R. has 
shown that for basic random processes represented by a continuous stochastic 
variate, every state between y and y + dy is entered only once—a rather 
surprising result, which on a closer examination is found to be satisfactory 
both mathematically and phenomenologically; F (y, t) is therefore given by 


FO.) = 5 O-t)+20+49}. (25) 


Since every state is entered only once, the probability 7 (y, 7) can be expressed 
neatly in terms of F (y, f), 
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a(y,) =f F (y, re (dr, (26) 


e~ (7) represents the probability of continuance in the state y. 


mn 1) = — An(y, t) + F(y, 0. (27) 


Substituting the value of F (y, t) as given by (25), we obtain the equation 
(17) thus proving the correctness of our expression for F (y, t). The problem 
of recurrence has no meaning in this process, since every state can be entered 
only once. 


t 
Generalisation of the Process | m (7) dr. 


We were concerned so far with an S.O.P. in which with each realised 
Poisson event, a value + | or — | was associated with equal probability 4. 
We can generalise this process by associating a value a with each event, a 
itself being random, with a probability frequency function ¢(a) which is 
symmetrical about a = 0. For such a process we ask for the probability 
distribution of the sum @ of the values associated with the Poisson points 
realised in an interval (0, t). It is easy to obtain the equation satisfied by 
the probability frequency function 7 (9, t) of this sum. 


(0,1) = f e*" Ade f$(a) (0 — a, t— 7)da +e 3(6), (28) 


To solve this we define p (u, t) and »(u) as the Fourier transforms of z (9, t) 
and ¢(a) respectively. Then 


p(u, t) = / e-a7 Adr p(u, t — 7) n(u) + eo, (29) 
A further Laplace transform with respect to t yields the L.T. of p(u, ¢) as 

p(t) = yz ON} + x15, (30) 
so that 

p(u, r) = I =a@) (31) 
Hence 


p(u, t) = exp [— A(1 — nw) ¢] (32) 
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and 


+00 


(0, t) = * f e¥6 exp [— A(1 — n(u)) t] du. (33) 


7 
—oo 


Thus the form of 7 (9, ¢) can be determined either analytically or numerically, 
t 
when the form of ¢ (a) is known. Again, denoting f 6dt by y, the probability 
0 


distribution of y can be found by the use of the regeneration point method: 


n(y,t)= [e** Adz [ $(a)a(y — a(t — 7); t— 7) da + ce 3(y) (34) 
and the Fourier transform 
p (u, t) of 7 (y, t) 


is given by 
p(u, t) = fer Adr ie (7) p(u, t — t) da + e~*. (35) 
Titbventiating the stilt equation, we obtain 
PCD) — — rp (u, 1) + Av (u, 1) 9 (ul (36) 
Hence 
p(u, t) = exp. D ft n (ut) — 1} dt] 
= exp. [- Mu) ; (37) 
where 
h(u) = / {1 — 9 (u)} du (38) 


The “‘ delta” approximation 


Though we have introduced the function ¢(a) we shall now show that 
the following approximation is possible and yields good results if ¢(a) is 
a function which rapidly falls off to zero as a differs from 0. To put it more 


+0 

precisely, assuming ¢(a) to bea symmetric function, f a*™41 4 (a) da =0 
+00 —co 

and / a2" 4 (a) da = « {a*"} exists, but we shall adopt the approximation that 


¢{a*} is very small, and its higher powers can be neglected. Since «{a?™} 
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is of the order of [¢ {a*}]" we shall assume all the moments are negligible 
except «{a*}. The same results can be obtained, if we assume 


$ (a) = ${ 8 (a — ap) + 5(a + ag)}, (39) 


where ay = 1/e {a®} and 5 is the Dirac-delta function. Hence in equation 
(34), if we adopt this approximation the generalised process reduces to the 
simpler type given by equation (16) which is explicitly solved. We shall 
call the above approximation the 5-approximation remembering that it is 
not ¢ (a) that it is approximated to a 5-function, which gives the trivial result 
that there is no dispersion at all, but that ¢ (ca) is a sum of two 5-functions as 
given in equation (39). 


4. APPLICATION TO MULTIPLE SCATTERING 


The above theory can be directly applied to the well-known multiple 
scattering problem. When a charged particle passes through matter, it 
suffers collisions with the: atoms of the matter. If we interpret the Poisson 
events to refer to these collisions and a to the angle of scattering in a single 


t 
collision then y(t) = / 6(7) dr is the lateral displacement of the particle 


provided we assume that 5’ a; = @ the total angular scattering, is small and 
é 


that the particle pursues almost a straight path with slight lateral displace- 
ment in the same plane (“ t — y”’ plane); tis the thickness of matter traversed. 
This problem has been elaborately investigated by Snyder and Scott and our 
equation (37) is identical with theirs. On the 5-approximation we get the 
explicit solution (22). 


5. MULTIPLE SCATTERING IN MULTIPLICATIVE PROCESSES 


The problem of multiple scattering of a single particle is simple when 
compared to the multiple scattering in multiplicative processes. Such a 
situation arises for example in a cascade, i.e., if we admit the possibility of 
new particles being produced on collision. The essential complication con- 
sists in the randomness of the points of birth of new particles and the direc- 
tions in which they are projected at the time of birth, since the lateral dis- 
placement of any particle is dependent upon the above factors. Many 
attempts have been made to deal with the three or two dimensional cascade 
spread. It is not difficult to write down formal integral equations using the 
regeneration point method. But these equations are intractable, and are 
not amenable to numerical computation. We here suggest three simple 
models whose stochastic features will be discussed in a later paper supported 
by numerical work, 
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Model 1.—We shall assume p dt is the probability that a particle collides 
and produces an-iher new particle and A dt is the probability that it collides 
and gets merely <cattered. When it produces a particle, the ‘ parent’ and 
the ‘ offspring’ continue to move in the same direction as that of the parent 
before collision. Defining P(n, y, tf) as the probability that there are n 
particles with displacement less than y, using the standard arguments of the 
regeneration point method, we can at once write the integral equation 


P(n, y, t) = fem = [A [P(n,y — Ot — 2, t— 6 (0) 8 
0 6 
+ p> uP (n,, y, t — 7) P (ns, y, t — 7)] dr 


nitna=n 
+ € OW TH (y) 8)", (40) 
where 6," = 0 if n #1 and 5,1= 1 and H(y) = 0 if y <Oand H(y) = 
if y >0. The mean number of particles with displacements less than y can 
be obtained easily as ' 


ext P(y, t) where P(y, 1) = / a (y’, t) dy’ 


and 7 (y, t) is the solution of equation (34). It must however be stated that 
the fluctuation about this mean number cannot be obtained so simply in 
terms of 7 (y, f). 


Model 2.—We shall assume the ‘ parent’ to move in the same direction 
and the ‘ offspring’ to be projected at an angle f to the direction of the 
‘parent’, where 8 has the probability frequency function %(f8) symmetric 
about zero. The corresponding integral equation of P(n, y, f) is given by 


P(n, y,t)= fe Om Adr[ [ P(n, y — A(t — 7), t — 7) $(8) dO 
é 8 


+ 5 [P(n,y— B(t— 7), t—7)P (ne, y, t— 7) ub (B) dB] 


nytm=n B 


+ ete)? A (y) 3,". " (41) 


Model 3.—The ‘ parent’ and ‘offspring’ are projected at angles f, 
and f, with respect to the original direction of the ‘ parent’ and f, and B, 
have the joint probability frequency function %(8,, 8.). The integral equa- 
tion for this model is given by 


Pin, y, t) = fe OM" Ade [ fP(n, y — O(1 — 2), 1 — 7)4(0) dd 
0 6 


+ 5 f[fP@,y—Bi(t— 7, t— 7) P (ms, y — Po(t — 2), t-7) 
mytnmg=n Br Ba 


< pib (By, Bz) dBydBe] 
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+ er Ott Hi (y) 8,", (42) 


We propose to discuss these equations and solutions under the delta 
approximation in a later contribution as stated before. 


Blatt (1954) has obtained similar equations, but he has used negative 
probabilities in his derivation. This, in the opinion of the authors, is not 
warranted, as no specific advantage is gained by the introduction of such an 


artificial device. A short discussion on this is included as an Appendix to 
this paper. 
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APPENDIX 


Let us consider equation (34). We have assumed that in an interval df, 
A dt is the total probability for scattering and A¢ (a)dadt, is the probability for 
scattering through an angle (a,a-+ da). 1 — Adt is the probability that no 
scattering occurs. Blatt introduces the quantity — Adt as the fictitious 
negative probability for forward scattering, i.e., through a =0O and says 
that the total cross-section for scattering is thus reduced to zero. The pro- 
bability of continuing without scattering which is 1 — Adt according to us, 
is taken to be 1 by Blatt. In other words, he replaces ¢ (a) by ¢(a) — 5 (a), 
5 being the usual Dirac delta function. (34) then takes the form 


a(y, t) = [ drdz [ {$ (a) — 3(@)} a (y — a(t — 7), t — 7) da + 8(y). 


Differentiating this expression with respect to ¢, we get exactly the same 
equation for 7(y, 7), as is obtained by differentiating (34). Hence, no 
simplification is obtained by splitting up the quantity 1 — Adt into 1 and 
— dAdt, and attributing — Adt to forward scattering and hence the principle 
of “‘ invariance” with respect to the addition of a delta function term with 
its singularity at a = 0 to 4(a) is without physical content. 
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DIETHOXY ANILINE 
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Received December 21, 1955 
The dipole moments of 2-amino-1:4-dimethoxy benzene and of 
2-amino-1 : 4-diethoxy benzene have been studied in benzene and in 
carbon tetrachloride solution. The values of the moments are discussed 
in relation to the moments of the unsubstituted parent carbon com- 
pounds. The interaction between three groups of similar characteristics 
is considered in relation to their steric effects. 


HYDROQUINONE and its methyl and ethyl ethers are well known instances® 
where the moments of two substituents in the para position do not cancel 
each other. It was felt that the introduction of a third substituent might 
show additional features similar to those of p-xylene derivatives. The pre- 
sent communication deals with the amino derivatives of the two quinol ethers. 


EXPERIMENTAL 


The solvents used were purified as in the previous communication on 
para-xylene derivatives.! 


2-aminol-1 : 4-diethoxy benzene: The Kodak sample was first decolour- 
ised by crystallisation from water in the presence of decolorising carbon 
and then recrystallised from distilled water. White flakes m.p. 86° C. 


2-amino-|: 4 dimethoxy benzene was also prepared in the same manner 
from a Kodak sample. White flakes, m.p. 81°C. 


The measurements of dielectric constant were by a ‘ heterodyne beat 
method’ as described earlier’ and refractive indices and densities were also 
measured similarly. Care was needed to exclude light and air during re- 
crystallisations of these compounds. 


The relevant experimental results are collected in Tables I and II below. 


DISCUSSION OF RESULTS 


In the parent compound, the existence of a dipole moment has been 
attributed to the presence of free rotation about the Cy,,j;—O bonds.? The 
introduction of a substituent group ortho to one of the ether groups can be 
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TABLE I 


Measurements with 2-amino-1:4-dimethoxy benzene 


at 32+0-01°C. 


























Weight Density Refractive Dielectric 
Solvent fraction of of index of constant of 
solute solution solution solution 
Benzene ase 0 0-86521 1 -4934 2°2588 
0:005759 0-90390 1-4942 2-2692 
0-01226 0-90478 1-4948 2°2831 
0-01766 0-90680 1-4956 2°2868 
0-02351 0-90738 1-4962 2:2978 
Carbon tetra- 0 1-5705 1-4536 2°2155 
chloride 
0-002941 1-6390 1-4540 2+2256 
0-004960 1-6376 1 -4547 2°2285 
0-006410 1 -6367 1-4550 2-2290 
0-007080 1-6362 1-4552 2-2300 
TABLE II 
Measurements with 2-amino-1:4-diethoxy benzene 
at 32+0-01°C. 
Solvent w d n € 
Benzene me 0 0-86521 1-4934 2-2588 
0-007344 0-90400 1-4935 Ss 
0-009795 0:90452 1-4942 2:2704 
0-01383 0-90500 14945 2°2835 
0-01783 0-99635 1-4948 2-2990 
Carbon tetra- 0 1-5705 1-4536 2°2155 
chloride 
0-002871 1-6348 1-4538 2°2236 
0-006354 1-6327 1°4541 2-2286 
0-007018 1-6323 1-4545 2+2338 
0-007533 1-6320 1-4548 2°2368 
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expected to modify this picture. Here, too, it is worth studying how far the 
different methods of calculation of dipole moments lead to the same value. 
The results are presented in Table III below. 
TABLE III , 
Dipole Moments of Dimethoxy aniline and Diethoxy aniline 





Halverstadt-Kumler 





Compound 





Pajcco Mx” oFe Moment 
2-Amino-1 : 4-dimethoxy 
benzene .. 67-89 53-35 14-54 0-85 
2-Amino-! : 4-diethoxy 
benzene eee fb 60-82 10-89 0-74 





It is not our present purpose to consider the relative merits of the different 
approximations used in the evaluation of dipole moments, but, most of the 
previous studies on related compounds have used the Halverstadt-Kumler 











method and for comparisons, these values alone will be used. Everard 
and Sutton*® have studied the dichloro and dibromo derivatives and these 
values together with those of the parent ethers are collected in Table IV. 
TABLE IV 

Compound Dipole moment Reference 
1:4-Dimethoxy benzene i - 1-81 (1-70) 4 (5) 
2:5-Dichloro-1: 4-dimethoxy benzene ee 3 (a) 
2: 5-Dibromo-! : 4-dimethoxy benzene «+, 3983 3 (b) 
2-Amino-| : 4-dimethoxy benzene s. 0°85 Present work 
1:4-diethoxy benzene... - ea 4 
2-Amino-1 : 4-diethoxy benzene -- 0°74 Present work 





The alkoxy groups are well known examples of groups having appre- 
ciable mesomeric moments and this would naturally impose a certain degree 
of restraint on free rotation. The presence of any substituent on the 2 and 5 


positions may be expected to further modify this and it is natural to expect 
Ad 
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TABLE V 
Moment 
Compound Moment in carbon 
in Benzene tetrachloride 
1: 4-Dimethoxy benzene a) - 1-70 1-81 
2-Amino-|: 4-dimethoxy benzene ba 0-85 1-37 
2-Amino-1: 4-diethoxy benzene .. = 0-74 1-21 








FIG 1 2-NH;1:4-C,H,OCH), 
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that the moments may be reduced by such substitution. This is clearly borne 
out by Table IV. The steric interaction which requires the groups to take 
up a multiplanar arrangement is also illustrated by Figures 1 and 2. 





@---+-5- 





FIG. 2. 24NH,1:4-C,4 (OCH), 
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In the parent ether itself, the presence of the two alkoxy groups in para 
position, by mutual action, may cause a reduction in the resonance effect and 
two extreme configurations are the cis and trans configurations. The former 
will be the more polar group but there is no compelling reason why the cis 
form should be favoured. Possibly, we are having here an instance of 
appreciable atom polarisation.® 


Several instances of interactions between ortho substituents of this type 
have been reported with the mono ethers.” From their study of the fluoro, 
chloro and bromo anisoles and of chlorophenetole, Anzilotti and co-workers 
find that the ortho substituents take up what they call a ‘ trans-trans ’ arrange- 
ment and a similar conclusion is reached by Smyth® from a study of di- 
tertiary butyl quinol dimethyl ethers (moment 1-42). In the presente instance 
also, a stable arrangemént requires that one of the groups moves in a direc- 
tion in the plane of the ring while the other moves in a direction perpendicular 
to the plane of the ring. 


While the unsubstituted ether has the preferential cis configuration of 
the methoxy or ethoxy groups, a trans configuration is favoured in the amino 
derivatives by the possibility of hydrogen bonding between the amino group 
and the oxygen of the methoxy or ethoxy group which is possible only in 
this: arrangement. Ortho methoxy phenols are known to indicate such 
interaction between the hydroxy group and the oxygen of the methoxy 
group.® 

CH, 
CH; 


18) 
Aww 


e 


cis form O 


CH; | 
\ ys OC,H; trans form 
CH, 


This spatial arrangement also results in a reduction in the effective 
moment of the compounds. The overall lowering is thus a cumulative one, 
the reduction in moment being greater in the ethoxy compound possibly by 
the stronger hydrogen bonding induced by the ethyl as compared to the 
methyl group. One consequence of the interaction is that the moment is 
reduced to a value below that of aniline. 
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ADDENDUM 


It will be noticed from the figures that the difference in moments between 
the solution in benzene and in carbon tetrachloride is very large. The 
abnormality is with the carbon tetrachloride value. Apart from the usual 
difference that has been generally observed with the change in dielectric 


constant of solvent, even though the solvent is nonpolar, fine needle-shaped 
colourless crystals which are not the same as the original compound dissolved 
have been consistently obtained. The nature and identity of the compound 
has not yet been established but the formation of this new compound is at 
least in part responsible for the abnormal values obtained. 
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INTRODUCTION 


ULTRASONIC velocities in pure liquids are known to be related to their densi- 
ties and compressibilities. Ultrasonic absorption is a more complex pheno- 
menon and several types of relaxation effects have been postulated as causing 
absorption besides such well-known agencies as viscosity and thermal con- 
ductivity. Even for pure liquids the experimental observations of ultrasonic 
propagation have not been fully and satisfactorily explained. Therefore 
when we pass on to liquid mixtures, the dependence of velocity and absorp- 
tion on concentration is not at all easy to predict. While some expressions 
have been proposed relating velocity with concentration of the mixture and 
with the properties of the constituents, no such attempt has been made in 
the case of absorption. Much experimental work has to be carried out 


before any generalisations on the ultrasonic absorption in liquid mixtures 


can be made. It is hoped that the results presented in this paper will have 
some value in this direction. 


The experimental method employed is the pulse technique of Pellam 
and Galt (1946) and Pinkerton (1949) with some modifications. Instead 
of a reflector, a second crystal has been used as receiver. A time mark 
generator and pulse generator both of which are controlled by a 100 Kc./s. 
crystal oscillator have been designed and constructed by the author for this 
purpose. As the method is by now well known no further description is 
given here. All observations were made at the Ist, 3rd, 5th and 7th harmonics 
of a pair of identical X-cut quartz crystals of fundamental frequency 
2-12 Mc/s. The temperature was controlled to within +0-5°C. The 
absorption measurements, which were confined to the Fresnel region of the 
ultrasonic beam, are made with an error not exceeding + 5%, and velocities 
with an error of + | m./sec. 


The liquids are all standard analytical reagents of guaranteed purity. 
RESULTS AND DISCUSSION 


The results in the five mixtures are given in Tables I to V. The velocity 
(V) is in metres/second and the absorption a’ = a/f? is in 10-™ db sec,? mr, 
106 
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Concentration is expressed as the percentage volume of one of the consti- 
tuents with respect to the total volume of the mixture. Some of the more 
interesting results are represented graphically in Fig. 1. 





(A) (8) 


Dioxane vol. Percent Ethyl Alcohol vol. Percent 
i i 
50 





(¢c) 


Fig 1: ULTRASONIC VELOCITY 
AND ABSORPTION IN 


(A) Elayé Methyl Ketone - Dioxane 


(8) Benzyl Alcohol ~ Ethyl Alcohol 


(Cc) Benzyl Alcohol - iso Propyt 
Alcohol. 


—B— VELOCITY 


—O— ABSORPTION 











& Propyl Alcohol vol. Percent 





79 
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Two of the five mixtures studied show normal behaviour in that their 
absorption and velocity move gradually from that of one liquid to that of 
the other. These two are benzene and ethylene dichloride, and ethyl alcohol 
and oleic acid. The ethyl methyl ketone and dioxane mixture shows a sharp 
and prominent peak, in its absorption curve and the benzyl alcohol-ethyl 
alcohol and benzyl alcohol—isopropyl alcohol mixtures show sharp dips 
in absorption curves. The results in the five mixtures will now be studied 
separately. 


(1) Benzene-Ethylene Dichloride—The velocity of the mixture decreases 
somewhat rapidly at the lower concentrations of ethylene dichloride in 
benzene but the greater part of the curve is a straight line. 

The absorption decreases very rapidly with addition of small quantities 
of ethylene dichloride to benzene as observed by Wada and Shimbo (1952) 
in the case of other benzene mixtures. This they ascribe to the more efficient 
molecular collisions in the presence of molecules of a liquid of low absorption. 


TABLE I 


Ethylene Dichloride and Benzene 
Temperature 33° C. 








Ethylene dichloride V a’ 
0 1272 829 

4-6 1358" 765 

6-25 1257 736 

9-1 1243 677 

16-7 1217 586 

25-0 1206 340 

50-0 1196 157 

66-7 1193 152 

75-0 1189 130 


83-37 1175 
100-0 1170 120 
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TABLE II 
Ethyl Alcohol and Oleic Acid 


Temperature 28° C. 











Cc 
Oleic Acid V a’ 

0 1150 49 
16-7 1209 106 
25-0 1220 165 
50-0 1279 225 
75-0 1344 412 
83-3 1357 a7 
87-5 1365 427 

100-0 1378 420 





(2) Ethyl Alcohol-Oleic Acid.—In this mixture oleic acid is a very 
viscous liquid and has a fairly high absorption coefficient. The velocity of 
the mixture as well as its absorption coefficient vary gradually from the 
oleic acid values to ethyl alcohol values. It is, however, interesting to 
observe that the effect of small quantities of alcohol even upto 25 per cent. 
by volume of the mixture is very little. The absorption is nearly constant 
in that region and the velocity gradient is much smaller than in the later 
portion of the curve. This can be ascribed to the high viscosity of the oleic 
acid which is not changed very much by the addition of small quantities of 
ethyl alcohol. 


(3) Ethyl Methyl Ketone-Dioxane.—This is a mixture of a highly polar 
liquid (the ketone) and a non-polar liquid. The former has also a low 
absorption compared with the latter. The same applies to velocity. 


The velocity of this liquid mixture varies linearly with concentration if 
the two points corresponding to peak absorption are ignored. Here there 
is a sudden break in the velocity-concentration curve and is shown by a 
dotted line. The behaviour of absorption also is similar. If we ignore the 
sharp peak, the curve is similar to the normal one for benzene-ethylene 
dichloride, That part of the curve between the foot of the peak and pure 
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TABLE III 


Ethyl Methyl Ketone and Dioxane 
Temperature 28° C. 








¢ 
Dioxane Vv a’ 

0 1165 43-0 
10-0 1182 33°1 
16+7 1193 60-3 
20-0 - 100-0 
25-0 1203 110-0 
26°4 1188 91-7 
37-5 1240 48-0 
50-0 1237 

57-9 1257 53-0 
75-0 1281 73-0 
83-3 1285 81-0 
100-0 1313 155-0 





dioxane represents the behaviour of a mixture of two unassociated liquids 
whereas the peak itself is similar to that in ethyl alcohol-water mixtures 
(Burton, 1948). In this latter case, the peak was ascribed to breaking up 
of the molecular association of water by alcohol. If we assume a certain 
amount of association for the ketone molecules in view of their high dipole 
moment and if, further, the dioxane molecules are assumed to break up this 
association, the peak can be easily explained. The curve for this mixture is 
dissimilar to that for water-ethyl alcohol only because here we are dealing 
with liquids of widely differing absorptions whereas for the latter mixture 
the absorptions of the pure liquids are not so different. 


(4) Benzyl Alcohol-Ethyl Alcohol, Benzyl Alcohol-Isopropyl Alcohol 
Mixtures.—In these two cases the velocities vary nearly linearly with con- 
centration but the absorption shows a sharp dip (reverse of a peak) at small 
concentrations of the aliphatic compounds in benzyl alcohol. This type 














of result has not been reported before although when unassociated liquids 
of nearly equal attenuation are mixed, they are known to pass through a 


gradual minimum of absorption. 
TABLE IV 


Benzyl Alcohol and Ethyl Alcohol 
Temperature 27° C. 
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Cc 
Ethyl Alcohol Vv a’ 
0 1508 82-3 
16°7 1457 124 
25:0 1436 66 
31-0 1422 140 
33-0 1396 140 
35-0 1386 140 
40-0 1383 115 
50-0 1351 107 
75:0 1245 74 
100-0 1156 49 
TABLE V 
Benzyl Alcohol and Isopropyl Alcohol 
Temperature 27° C. 
¥ Cc 
Propyl Alcohol Vv a’ 
0 1508 82-3 


4-0 
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Both the constituents of the mixture are polar molecules having low 
absorption. Following Hall (1948), if we ascribe the excess absorption of 
the pure liquids to a structural relaxation of the associated molecules, it means 
that so far as these frequencies are concerned the equilibrium between 
external and internal degrees of freedom of single molecules adjusts itself 
very rapidly and hence does not contribute to absorption. The dip in the 
absorption of the mixture can be explained then by assuming that there is 
a minimum association at that point (or perhaps none at all) and the indivi- 
dual molecules propagate ultrasonics without much absorption. On either 
side of the dip, molecules again associate giving rise to the usual absorption. 


SUMMARY 


Ultrasonic velocities and absorption in five binary liquid mixtures have 
been studied in relation to their concentration. They are benzene-ethylene 
dichloride, ethyl alcohol-oleic acid, dioxane-ethyl methyl ketone, benzyl 
alcohol-ethyl alcohol and benzyl alcohol-isopropyl alcohol. The rapid 
increase in absorption of benzene with the addition of small quantities of 
ethylene dichloride, the negligible influence of small quantities of ethyl alcohol 
on the high absorption of oleic acid, the sharp peak in absorption-concen- 
tration curve in the case of ethyl methyl ketone-dioxane mixtures, the sharp 
minima in the last two mixtures are explained qualitatively. 


In conclusion, the author wishes to thank Professor S. Bhagavantam 
for his guidance and keen interest in this work. 
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INTRODUCTION 


THE measurement of the vertical intensity of cosmic radiation has been 
carried out by different workers—® upto a maximum depth of about 
3,000 metres water equivalent (mwe). Barrett et al.* have combined the 
various measurements and after suitable normalisation, have found it possible 
to draw a smooth curve passing through all the experimental points. 


Recently, we have made measurements of vertical intensities at depths 
100, 381, 475, 684 and 885 mwe** in the Kolar Gold Mines (Mysore State). 
Our points also fall very close to the composite intensity-depth curve of 
Barrett et al. (Fig. 1). The slope of this composite curve is found to increase 
continuously with depth. It does not exhibit any abrupt changes in the 
slope, as was believed to be the case, a few years ago. 


Absorption measurements of Barrett et al. at a depth of 1,574 mwe, 
have clearly shown that the penetrating particles observed underground, are 
not locally produced in the rock above the recording apparatus, but are 
produced high up in the atmosphere, and traverse all the way through the 
tock before reaching various depths. Cloud chamber experiments, carried 
out below ground, have revealed that the general features of the penetrating 
particles below ground are very similar to those exhibited by the penetrating 
particles at sea-level. This, together with the fact that the least interacting 
particle known so far, is the u-meson, makes it quite certain that the under- 
ground penetrating particles are essentially the high energy »-mesons en- 
countered at sea-level. 


Therefore, the behaviour of the intensity-depth curve depends upon 
(a) the energy spectrum of y-mesons at sea-level, and 
(b) the range-energy relation of u-mesons. 





* The depths quoted in (3) were in crror and have been corrected in this paper. 
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It is well known that p-mesons are not produced directly in nuclear 
interactions. They are the decay products of heavier mesons. Therefore, 
the sea-level spectrum of u-mesons depends upon 
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(i) the production spectrum of the parents of u-mesons, and 


(ii) the mass, mean life-time, absorption mean free path, and relative 
abundance of the parent particles. 


At present two definite sources of u-mesons have been established, viz., 
the 7-mesons and some of the K-mesons. The mass, life-time and absorp- 
tion mean free path of z-mesons are very well established. The properties 
of K-mesons are however not so well determined yet. Among the K-mesons, 
which decay into ~-mesons, two types have been recognised. 


(a) the K,,, which undergoes decay according to the two-body decay 
scheme 


K,, > + v and 


(b) the K,,, or the kappa-meson which decays according to the three- 
body decay scheme 


K,, >a +?+2 


Evidence for K,,, undergoing a two-body decay, was first obtained by 
the Ecole Polytechnique group.* Bridge et al.5* measured the range of 
yemesons emitted by K,, mesons stopped in a multi-plate cloud chamber, 
and deduced the mass of K,, as 950+ 15 me. From direct measurements 
on K,, particle, as well as measurement of the range of the secondary 
pemeson, Armenteros ef al.5© obtained 935 + 15 me as the best mass value 
of Ky. Direct timing measurements on the life-times of K-particles by 
Robinson,* and Barker et al.’ indicate that the K,, life-time is of the order 
of 10-® secs. 


Evidence for K,, or kappa mesons was first found by Menon and 
O’Ceallaigh® in photographic emulsions. Direct evidence for the three- 
body decay scheme involving a »-meson is obtained from four cases (Bristol, 
Rome, Rochester and Paris), in which the secondary energies suggested a 
continuous energy distribution. The secondaries were identified as u-mesons 
by p-e decay. The best estimate of the mass of K,, is 1035 + 25 me.® 
Armenteros et al.’ give a life-time of the order of 10~® secs. for K,, particles. 
The relative abundance of K,, and K,, is not well known. An absolute 
upper limit to the ratio of charged heavy mesons to charged 7-mesons, pro- 
duced at very high energies, is deduced to be 25%.4 

At low energies all the 7 and K,-mesons decay into y-mesons. But as 
the energy increases, due to relativistic increase .of the life-time competition 
between absorption and decay sets in. This results in a progressive defi- 
ciency in the number of u-mesons at the high energy end. Using the decay 
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and absorption constants, it is possible to work out theoretically, the sea-level 


#-meson spectrum provided the production spectrum of the parent particles 
is known. 


In Section I of this paper, we have presented arguments to show that 
a measurement of zenithal angular distribution of penetrating particles at 
a depth of about 100 mwe, will enable us to accurately determine the pro- 
duction spectrum of the parents of u-mesons. We have also shown that 
this determination is fairly independent of the nature and composition of 
the parent particles. 


In Section II, we have given the results of an experiment on the angular 
distribution of penetrating particles, carried out at a depth of 103 mwe, 
with a hodoscoped counter telescope. The experimental results indicate a 
power law production for the parents of the ~-mesons with an exponent 
— 1-92 + 0-08. 


In Section III, we have calculated the theoretical intensity-depth curves 
for different conditions of production of the parents of u-mesons and also 
using different range-energy relations of y-mesons. The theoretical curves 
have been compared with the experimental intensity depth curve. The 
conclusions that can be drawn from the comparison are discussed. 


SECTION I. RELATION BETWEEN THE ANGULAR DISTRIBUTION OF 
UNDERGROUND u-MESONS AND THE PRODUCTION SPECTRUM 
OF PARENTS OF u-MESONS 


Let us consider, for the present, the production of 7-mesons alone in 
the nuclear interactions of the primaries. The 7-mesons are known to decay 
with a life-time 2-56 x 10-* secs. into u-mesons. There is now substantial 
evidence to show that the 7-mesons'* and primary protons which produce 
them interact with nuclear geometric cross-section. upto energies of the order 
of 1044+ev. Considering the competition between decay and absorption of 
m-mesons, it is easy to derive (see Barrett et al.) the following expression 
for the sea-level intensity of ~-mesons arriving at an angle @ with respect to 
the vertical, with an energy greater than Eg at production: 


1(Es, ) =K, - (22) Be 
o” (€+ 1) E,cos 0+ «B, + 





VBE, cos 0 (1) 
ze +1 
where E, = minimum energy of the u-mesons at the production layer. 


« = the exponent in the integral spectrum of 7-mesons at pro- 
duction, assumed to follow a power law. 
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“ m, c? .¥,.H 
CT 
m,,c? = rest energy of 7-mesons = 0-140 Bev. 

Yn = energy degradation factor in 7-u decay = 0-76. 

Tz = mean life-time of 7-mesons = 2-56 x 10-* secs. 

H = equivalent atmospheric depth in cm. = 6-46 x 10° cm. 

K, =a constant. 
Upto fairly high energies, about 100 Bev, it is reasonable to assume the range- 
energy relation of relativistic .~mesons to be linear, so that 

E=a.R (2) 

where R is the range of u-mesons in mwe, and a is a constant giving the rate 
of loss of energy per mwe. 


The minimum energy E, that a u-meson should have at the top of the 
atmosphere, in order to arrive at a depth D, at an angle @ with respect to 
the vertical, 


E,=aR, (3) 
where R, is the amount of matter encountered by the u-meson in the direc- 
tion @ and is given by 

Rg = Dicos 8. (4) 
From (3) and (4), 

Eg = a D/cos @. (5) 
Substituting the expression for Eg given by (5) in (1), we obtain the expression 


for the intensity of ~-mesons arriving at an angle @ with respect to the vertical, 
at a depth D mwe below the top of the atmosphere, 


Ip (9) = K,, - (SS 2) = 
(e+ 1)oD + «B+ 








VED (6) 
2e+1 


= constant. cos€ 6 (7) 


since y,, T;, D, B, and e are all constants. Therefore, the exponent in the 
cosine law for the angular distribution of underground yu-mesons gives the 
exponent in the generation spectrum of z-mesons. 


Let us now take into consideration the production of K,-mesons as well. 


If the production spectrum of K,’s and 7’s follow the same power law, then 
from an exactly similar treatment as above, we have 


Ab 
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Ip (0) =K (ese Br 





ys .D 





2 tie ("eS cos 6\€ Br 
(e+ 1)aD + «Be + VBeaD 
2e+ 1 
= constant. cos€ 0. 


where Kx, By, x, yx and 7, are constants referring to K,-mesons. 


Therefore, provided the production spectrum of K,-mesons and 7-mesons 
is the same, the angular distribution experiment gives the exponent in the 
production spectrum. 


In the above analysis, we have assumed the range-energy relation of 
p-mesons to be linear. So it is essential that the angular distribution experi- 
ment is carried out at a depth for which this is true. The range-energy 
relation given by George*® taking into consideration the energy losses due to 
ionisation, bremsstrahlung, pair production, nuclear interactions, and 
emission of Cerenkov radiation, is non-linear and can be expressed as follows: 


R’ = 1/b in(i + bE/d) (8) 
where R’ is the range of u-mesons expressed in gm. cm.~? and E is the energy 


of the 4-mesons expressed in Mev. b = 5-26 x 10-6 for rock of mean atomic 
number 10-5 and mean atomic weight 21. d= 2-65 Mev. gm-!cm-? 


In order to see to what extent the non-linearity of the equation (8) 
affects the validity of equation (7), at a depth of 100 mwe, values of E for 0 
varying from 0° to 75° were calculated from (8) and substituted in (7). It 
was found that the value of the constant in equation (7) varied only by 4% 
in going from the vertical direction to an angle of 75°. 


Therefore the measurements: of the angular distribution of penetrating 
particles should give a fairly accurate determination of the production spec- 
trum. Another advantage of working at a depth of 100 mwe is that the 
effect of the u-e decay in the atmospheres of these energies is negligible. 


in the following section we have described an experiment on the angular 
distribution of penetrating particles, carried out at a depth of 103 mwe. 
SECTION II. ANGULAR DISTRIBUTION OF PENETRATING PARTICLES AT 
A DEPTH OF 103 MWE 


The experiment on the angular distribution of penetrating particles 
was carried out at a depth of 103 mwe, in the Gold Mines at Kolar, where 
the vertical intensity of cosmic rays was previously measured upto a depth 
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of 885 mwe. The counter geometry employed for the experiment is shown 
in Fig. 2. The counters were used in each of the trays A, B, Cand D. Ten 
centimetres of lead were placed between the trays B and C. The separation 
between the extreme trays was 42cm. The counter telescope was mounted 
in a steel frame which could be tilted about a horizontal axis and set at any 
desired inclination with respect to the vertical. 
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Fic. 2. Counter telescope used for the angular distribution exp2:imt. 

With a hodoscope arrangement, similar to the one described in a pre- 
vious paper (reference 3), the following counting rates were recorded for 
inclinations of 0°, 45°, 60°, 65°, 75° and 90° with respect to the vertical. 

(a) two-fold coincidence rates A; D; of the nine parallel 
single counter telescopes formed by correspond- 


ing pair of counters in trays A and D (i= 1......9) 

(5) three-fold coincidence rates A; BD; | ed een 9) 

(c) three-fold coincidence rates A; CD; | ee pores 9) 

and (d) four-fold coincidence rates A; BCD; 3 eee 9) 


In Table I we have given the coincidence rates (d) and the sum of the coinci- 
dence rates (b), (c) and (d) in all the nine channels, for different inclinations 


of the telescope. These counting rates are corrected for showers as described 
below. 


Shower Correction 


As is the usual practice in all angular distribution measurements, the 
shower correction has been obtained by recording the counting rate with 
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the telescope tilted horizontal. For large angles, for which the shower 
correction is appreciable, the maximum contribution to the showers will 
be from the vertical direction. So far as these vertical showers are con- 
cerned, the disposition of the counters will remain practically the same for 
all angles greater than 45°. Therefore, the horizontal counting rate should 
give a fairly accurate shower correction, for all large angles. To verify 
experimentally whether the horizontal reading gives a fairly reliable shower 
correction, measurements were made at sea-level, where the effect of side 
showers on telescope counting rates is not as significant as underground. 
With the same telescope at sea-level, the ratio of four-fold to three-fold 
coincidences before, and after shower correction, were almost equal (0-71), 
whereas below ground, the two ratios were significantly different (0-65 and 
0-70 respectively), thereby showing that the contribution of side showers 
underground is more than at sea-level. The ratio of the number of particles 
causing four-fold and three-fold coincidences is a geometrical factor, arising 
purely out of leakage spaces in counter trays B and C. This factor should 
be the same at sea-level and underground. The ratio of four-fold to three- 
fold coincidences obtained after subtracting the horizontal reading, is found 
to agree with the geometrical factor at sea-level. This comparison indicates 
that the shower correction applied here is quite reliable. 


The counting rates at different inclinations, corrected for showers by 


subtracting the horizontal counting rate from the observed counting rates, 
are given in columns 5 and 8. The sum of the counting rates (b), (c) and (d) 
corrected for showers is plotted against cos @ on a log-log scale in Fig. 3. 
It is seen that a straight line can be fitted to the experimental points. The 
slope of the line fitted by the method of least squares comes out to be 
1-92 +0-08. Therefore, the angular distribution at a depth of 103 mwe 
can be expressed by the cosine law, 


I (8) = 1 (0) cos?-92+0-%9 
or the exponent in the production. spectrum of the parents of u-mesons is 
— 1-92+ 0-08. 
SECTION III. COMPARISON BETWEEN THE THEORETICAL AND 
EXPERIMENTAL INTENSITY-DEPTH CURVES 


Using equation (1) in Section I, and the production spectrum E--%2, 
we have worked out the theoretical intensity-depth curves, considering the 
production of 7-mesons alone in nuclear interactions of primaries, and 
assuming | 
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Fic. 3. Angular distribution of penetrating particles at 103 mwe. 


(a) the range-energy relation of u-mesons to be linear upto 300 Bev, 
and a constant energy loss of 2-2 Mev gm. cm.? 


(b) the range-energy relation given by George (equation 8 in Section 
Il), which takes into account the various processes by which 
#-mesons lose energy. 


These theoretical intensity-depth curves are compared with the experi- 
mental intensity-depth curve in Figs. 4 and 5, after normalising the theoretical 
intensity with the experimental value at a depth of 100 mwe. It is seen that 
while the theoretical curve drawn with the linear range-energy relation agrees 
reasonably well with the observed curve, the one drawn using the range- 
energy relation of George, deviates considerably from the experimental 
curve, The discrepancy is not removed to a great extent even if we take into 
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Fic. 4. Comparison between theoretical and experimental intensity-depth curves. 


(1) Experimental intensity-depth curve. 


Theoretical intensity-depth curves 
with e = 1-92. Linear range-energy relation (energy loss = 2:2 Mev. gm.~) cm.2 
and assuming 


(2) Only x’s are the parents of u-mesons. 
(3) Parents of y’s are ’s and Ky’s in proportion 3:1, m,u = 976 m., Tu = 10-* segs, 
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Fic. 5. Comparison between theoretical and experimental intensity-depth curves. 
(1) Exp2rimental intensity-depth curve. 


Theoretical intensity-depth curves 


with «= 1-92 and assuming only x’s are parents of u-mesons, and George’s 


range-energy relation (eq. 8). 


(2) Including nuclear losses"(b = 5-26 x 10-8 in eq. 8). 


(3) Excluding nuclear losses (6 = 2-80 x 10~ in eq. 8), 
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Fic. 6. Compariscn between theoretical and experimental intensity-depth curves. 


(1) Experimental intensity-depth curve. 
Theoretical intensity-depth curves 
with e = 1-92 and assuming George's range-energy relation (eq. 8) and ’s and 
Ky’s are parents of y’s in proportion. 
(2) : Ku = 3:1 (tye = 10° secs.). 
(3) a: Ku =1:1 (tu = 10-° secs.), 
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consideration only the ionisation and electromagnetic energy losses of 
p-mesons, and ignore completely energy loss due to nuclear interactions, 
as shown in Fig. 5 (curve 3). In order to see whether the production, in 
addition to 7’s of mesons other than 7’s which have a life-time less than that 
of a’s, as for example the K,’s, will remove the discrepancy between the 
observed and theoretical intensity-depth curves, we have calculated the 
theoretical intensity-depth curves for the following composition of the parents 
of u-mesons: 


(i) 7-mesons: K,-mesons = 3: | Tk, = 10° secs. 
(ii) 7-mesons: K,-mesons = 1: 1 Tk» = 10° secs. 
assuming rz, the mean fraction of energy of K, transferred to u-meson as 0-5, 


It is seen from Fig. 6, curve 2 that even a high percentage of K, produc- 
tion assumed in (i) is inadequate to explain the intensity-depth curve. Curve 
3 in Fig. 6 indicates that unless K,-mesons are produced in at least as large 
number as 7’s and have a life-time less than 10~° secs., the observed intensity- 
depth curve cannot be explained, if we take the integral production spectrum 
E-!-82 and take into consideration all the energy losses of u-mesons in the 
formulation of the range-energy relation of u-mesons. 


DISCUSSION 


If we consider the production of z-mesons alone in the nuclear inter- 
actions of the primaries, then in order to remove the discrepancy between the 
observed and theoretical intensity-depth curves, either 


(a) the range-energy relation of p-mesons has to be linear upto 
300 Bev. (energy loss not greater than 2:2 Mev gm.~! cm.?). 


or (b) the production spectrum of 7-mesons should be a power law with 
an exponent very much lower than 1-8 (about 1-5). 


There is sufficient experimental evidence to indicate that the energy loss 
of w-mesons due to electromagnetic processes like radiation and pair pro- 
duction, does increase with increasing energy of u-mesons. The increase of 
the soft component in equilibrium with the penetrating component is a clear 
indication of increased energy losses. Therefore there seems to be no 


reason for questioning the validity of the range-energy relation of u-mesons 
given by George. 


On the other hand, none of the angular distribution measurements 
carried out so far has given an exponent lower than 1-8 for the production 
spectrum. Follet and Crawshaw working at a depth of 60 mwe, obtained 
for the power of the cosine law the value 2:02 + 0-2, Quercia and Rispo]j} 
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have made a measurement of the angular distribution of penetrating particles 
at a depth of 40 metres in water. When we calculate the power of the cosine 
law from the intensities given by them for the various directions, we obtain 
the value 3-4. Again the sea-level measurements of Caro, Parry and 
Rathgeber!® on the energy spectrum of u-mesons, gives for the differential 
spectrum an exponent 3-0+ 0-2 for energies greater than 20 Bev. This 
again would be consistent with an integral production spectrum of 7’s with 
an exponent 1-8. 


If we consider, in addition to 7’s, the production of heavier short-lived 
mesons, then as pointed out already, the heavy mesons should be produced 
in as large a number as 7’s and should directly decay into u-mesons with 
a life-time less than 10-° secs., if the discrepancy between the observed and 
calculated intensity-depth curves is to be removed. However, the existing 
results on the production of heavy mesons. show that the ratio of production 
of ’s to Ky’s is greater than 3, and tiat the life-time of Ky’s is of the order 
of 10-* secs. The other possible explanation is that 7-mesons do not inter- 
act with the same cross-section as protons. But recent experiments with 
photographic emulsions have shown that 7-mesons do interact with a nuclear 
geometric cross-section upto 10" ev. There is no reason why the cross- 
section for the interaction should suddenly drop at higher energies. 


Calculation shows that if about 5% of the mesons produced in nuclear 
interactions are p-mesons, or short-lived particles decaying into u-mesons, 
then the discrepancy mentioned above is removed to an appreciable extent. 
However, there is no conclusive experimental evidence either in favour or 
against such a hypothesis, at present. 


SUMMARY 


The production spectrum cf the parents of the underground penetrating 
particles has been determined from a measurement of the angular distri- 
bution of the penetrating particles at a depth of 103 mwe. The results indi- 
cate a power-law production spectrum with an exponent — 1-92 + 0-08. 


This production spectrum has been used to deduce the theoretical in- 
tensity-depth curves for different compositions of the parents of u-mesons. 
A comparison of the theoretical curves with the composite intensity-depth 
curve leads to the following conclusions: 


(a) The experimentally determined intensity-depth curve cannot be 
quantitatively explained, if 7-mesons produced in the nuclear interactions 
of the primaries, are the only parents of u-mesons, and if they haVe the follow- 
ng properties ; 
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(i) the production spectrum is a power law with an exponent greater 
than 1-8. 


(ii) 7-mesons decay into »-mesons with a life-time 2-56 x 10- secs, 
(iii) 7-mesons interact with the same cross-section as the primaries. 


(6) The production, in addition to 7’s of other heavy mesons, which decay 
into u-mesons, will not remove the discrepancy between the observed and 
theoretical intensity-depth curves, unless these heavy mesons are produced 


in as large numbers as the 7’s and have a life-time for -decay less than 10-° 
secs. 


Existing experimental evidence indicates a ratio not smaller than 3 for 
the ratio of production of 7’s to heavy mesons and a life-time of the order 
of 10-* secs. for the decay of heavy mesons into p-mesons. 


Calculation shows that the discrepancy is removed to a considerable 
extent if about 5% of the mesons produced in nuclear interactions are 
p-mesons or short-lived particles decaying into w-mesons. At present there 
is no experimental evidence either in favour or against such a hypothesis. 
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EXCITED STATES OF C” IN THE REGION 11 TO 
14 MEV. FROM B" (d, 2) C” REACTION 


By Evanit KONDAIAH 


(Tata Institute of Fundamental Research, Bombay) 


Received January 13, 1956 


(Communicated by Prof. B. Peters, F.A.sc.) 


JOHNSON”? reported a number of excited levels in C!* in the region 11 to 14 
Mev. However, he did not use a separated B™ target. In the present study, 
a separated B™ target has been bombarded with 1 Mev. deuterons produced 
by the Cockcroft-Walton accelerator at the Australian National University 
in Canberra.* 200 micron thick C, nuclear emulsion plates have been 
exposed to the outcoming neutrons at 0° and 90° to the deuteron beam 
direction. Similar exposure with a blank target has also been made. While 
scanning the plates, only such recoil proton tracks that made an angle not 
greater than 10° with the incident neutron direction have been chosen; hence, 
the errors caused by uncertainties in dip measurements are less than 2%. 


The data obtained has been statistically analysed as follows: 
A polynomial regression of fifth degree, namely, 
Y =a-+ bX + cX? + dX® + eX*-+ fX5 


has been fitted to the observed data.* In this equation, Y is the number of 
tracks to be expected at the energy X, if fifth degree polynomial were to fit 
the entire data. In most of the cases considered, the fifth degree polynomial 
is found to take into account mainly the basic downward trend of the 
observations. If a peak is present in the expected values, it is taken as indi- 
cative of genuine maximum provided the observed value at that energy is 
above this peak. Where there is no peak in the expected values, X? test is 
used to locate the peaks in the observed data. Applying the X? test, we 
find that some of the values in the observed data are decidedly higher than 
the corresponding expected values and have highly significant X* values. 
If these few significantly high X* values are struck off from the total X*, the 
rest of the regression fits well by X* test. The conclusion, therefore, is that 
these few high values in the observed data represent the only energy ranges 
where maxima are seen to lie as far as the present data is concerned. 


* Refer to ‘Statistical Tables for Biological, Agricultural and Medical Research,” by Fisher 
and Yates. 
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Figs. | and 2 show the neutron spectra at 0° and 90° respectively. The 
positions of the neutron peaks located by the above-mentioned statistical 
analysis are shown with arrows. In both the figures we see the neutron 
continuum due to three and four particle break up of the compound nucleus 
C18, The neutron peaks corresponding to the different excited levels in C! 
are seen superposed over this neutron continuum. The shaded portion at 
the bottom of the two figures shows the background. 


TABLE I 


| 
(E,)o° | (E,)ooe Qooe | Qaverase 


! 


0-85 +0-05 | —0-13 | 0-6 0-08 | —0-18 | —0-15 13-9 +0°1 








1:05+0-05| 0:06 | 0-85+0:05 | 0-09 0-08 | 13-6+0+1 


| | 
1-95 + 0-05 0-95 a | ie a 12-9 +01 








2°75 + 0°05 | 1-77 2°25 + 0-05 re | 1-69 12:0 + 0:1 





Table I summarises the results obtained from this investigation. The 
Q values are calculated from 


13 1 — 5 
Qo = 75 En—- & V 2En Ep - cos 6—<€ Ep, 
where 9 = 0° or 90° as the case may be; Ey, is the energy of the neutron peak 


and Ep is that of the bombarding deuteron beam. Q,,,.,4 is taken as 13-724 
Mev.” ‘ 


The neutron peak corresponding to the 1-95 Mev. peak in the 0° spec- 
trum should be situated at about 1-65 Mev. in the 90° spectrum. Though 
there is some rise in this region of the 90° spectrum, this does not appear 
to be significant in the statistical analysis. 


I wish to express my thanks to Professors D. D. Kosambi and B. Peters 
for discussion. Mr. U. V. Ramamohan Rao’s help in carrying out the 
statistical analysis is gratefully acknowledged. Part of this work was done 
during my stay at Canberra and I am thankful to Prof. E. W. Titterton for 
his kind co-operation. 
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POLAROGRAPHIC BEHAVIOUR OF METALS 
IN ETHANOLAMINES : 


Part I. Iron 


By R. S. SUBRAHMANYA 


(Department of General Chemistry, Indian Institute of Science, Bangalore-3) 
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(Communicated by Prof. K. R. Krishnaswami, rF.A.sc.) 
INTRODUCTION 


WoLFSON! AND Jessop” have mentioned that iron gives a polarographic wave 
in alkaline triethanolamine solutions, but no detailed investigations have 
been reported. In the present work mono-, di-, and triethanolamines have 
been employed as cationic complexing agents to study the polarographic 
behaviour of the ferric iron. 


EXPERIMENTAL 


Monoethanolamine (B.D.H. reagent grade) and pure di- and triethanol- 
amine (Fisher Scientific Co.) are employed in the base solution. Acetic 


acid and sodium hydroxide solutions have been used to regulate the pH of 
the base solutions. No maximum suppressor is used in base solutions con- 
taining triethanolamine but in solutions containing mono- and diethanol- 
amines 0-005% gelatin is used. 


Details of procedure for running and evaluating the polarograms have 
been described: previously.** The values of ‘m’ for the. dropping . mercury 
electrodes used are 1-370, 1-227, 1-120 and 1-150mg. per second, while 
the droptime has been adjusted to 5 seconds a drop in 1N potassium 
chloride. 

RESULTS 
(A) Base Solutions Containing Triethanolamine 
1. Effect of pH. 


In these experiments acetic acid and sodium hydroxide are employed 
to adjust the pH of the solutions. In solutions of pH less than 11, the con- 
centration of sodium ion is maintained at 0-2 M. ,The colour of the solutions 
below pH 9 is yellow, while above this value it is colourless. The results 
are given in Table I from which it is clear that above pH 9-8, the polarograms 
exhibit two waves at about — 0-9 V and — 1-5 V vs. S.C.E. corresponding 
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to the reduction of iron (IID) to iron (II), and iron (I) to metallic iron res- 
pectively. Below pH 9-8, the second wave has a tendency to get merged 
with the final current rise. The height of the second wave is double that of 
the first wave between pH 9-8 and 11, and not at higher values. The half- 
wave potential becomes less negative with a diminution in pH. The values 
of E;—E, indicate that the second wave is irreversible at all pH values, while 
the first wave is irreversible only between pH 8 and 11. 


TABLE I 
Effect of pH on the polarographic behaviour of iron in triethanolamine* 


Conc. of iron (III) 1-182 millimolar; Conc. of triethanolamine 0-10 molar 








(Ej) (E,)s Heights of the waves 
pH VS. vs. (E;—E,) (Ez—E,)2 in terms of divisions 
S.C.E. S.C.E. of the graph (S/70) 
(Volt) (Volt) (Volt) (Volt) Ist wave 2nd wave 
13-50 —1-072 —1-675 0-060 0-050 11-0 9-8 
12-70 —1-034 —1-+593 0-060 0-060 12-6 15-9 
11-00 —0-870 —1-470 0-083 0-085 12-0 24-0 
9-80 —0-695 —1-370 0-085 0-080 12-0 24-0 
8-30 —0:422 Nowave 0-086 No wave 11-8 No wave 
8-05 —0-396 e 0-089 os 12-0 i 
7-10 —0-363 a 0-061 eh 7:3 ~ 
6°85 —0-316 eB 0-057 i 6:9 s 
6°40 —0-256 o 0-054 * 7:3 me 
5-80 —0-189 e 0-045 - 8-0 = 
3-70 Very close to 
+0-30 





* In this table E4 represents the half-wave potential and the subscripts 1 and 2 refer to the 
first and the second waves respectively. 

It has been noticed that interesting data regarding complex formation 
and polarographic waves of good analytical value are obtained in the alkaline 
region. Hence a detailed study employing sodium hydroxide, sodium 
carbonate, acetate, and ammonium chloride—ammonium hydroxide has 
been made and the results are given below. 


2. Effect of sodium hydroxide, sodium carbonate and ammonium chloride- 
ammonium hydroxide on the ratio of the heights of the two steps, the 
half-wave potential and the reversibility of polarographic waves in pre- 
sence of triethanolamine. 


(a) Sodium hydroxide-—Some of the polarograms showing the effect of 
sodium hydroxide at different concentrations of triethanolamine are given 












a i ae ele ee eee ae ee)” ae ee, 
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in Fig. 1. The average value of the results of analysis of the waves are given 
in Table II. The results indicate that for a given concentration of tri- 
ethanolamine, the half-wave potential of the first wave increases in the begin- 
ning and then tends to remain constant, as the concentration of sodium 
hydroxide is increased. When the concentration of sodium hydroxide is 
kept constant (below 1-0 molar) the half-wave potential of the first wave 
practically remains constant for low concentrations of triethanolamine 
(below 0-1 molar) and becomes more positive with an increase in the con- 
centration of sodium hydroxide. In a good number of cases, these values 
could not be determined on account of the interference of the final steep 
current rise. It is interesting to note that in most cases the values of 
(ig)e/(ia), is less than 2. 


TABLE II 


Effect of sodium hydroxide on the polarographic waves of iron (IIT) 
in presence of triethanolamine 


Conc. of iron (IIJ) 1-182 millimolar 








Conc. of Conc. of =-Half-wave Half-wave Slope of 

sodium triethanol- potential potential log plot* (ig)o/(ia)s 

hydroxide amine Ist wave 2nd wave 2nd wave 

(E;) vs. S.C.E. 

(Molar) (Molar) (Volt) (Volt) (Volt) 
0-10 0-002 —1-029 —1-530 0-070 1-56 
0-10 0-020 —1-030 —1-616 Tt 1-31 
0-10 0-100 —1-034 —1-+593 Tt 1-21 

‘0-10 0-300 —1-018 wie 
0-10 0-500 —1-013 
0-10 0-700 —1-004 
0-10 1-000 —1:00 ... ; 
1-00 0-002 —1-069 —1-579 0-035 1-81 
1-00 0-020 —1-074 —1-641 0-045 1-42 
1-00 0-100 —1-072 —1-675 0-045 0-83 
1-00 1-000 —1-058 —1-733 Tt 0-23 
5-00 0-002 —1-072 —1-414 0-045 1-77 
5-00 0-020 —1-079 —1-450 0-038 1-74 
5-00 1-000 —1-076 ‘ie 





* Slope of the log plot of the first wave is about 0-060 Volt. Some log plots are given 
in Fig. 2. 
+ Waves not well formed. 


(b) Sodium carbonate.—The results are summarized in Table III. The 
two reduction waves are very well developed and the ratio of the heights is 
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1: 2 in contrast to the behaviour noticed in. base solutions containing sodium 
hydroxide.. The half-wave potential shifts to less negative values with an 


0-60 


0-40}. 


0:20 | 


0-00}- / 2 
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Curves 1 and2 were obtained in base solutions containing 1M triethanolamine + 0-1M 
sodium hydroxide and 0-02M triethanolamine + 0°1M sodium hydroxide respectively. 
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increase in the concentration of iron. At a constant concentration of ferric 
iron, the half-wave potentials of the first and the second waves progressively 
increase with an increase in the concentration of triethanolamine. The 
values of E;—E, indicate that at low concentrations of triethanolamine 
(0-02 M), the first wave is reversible while at higher concentrations of tri- 
ethanolamine it is irreversible. The second wave, however, is practically 
irreversible. It is surprising that a reversible value (0-029 V) should have 
been obtained with 0-1 M triethanolamine and 0-1 M sodium carbonate, 
when the concentration of iron is 1-5 millimolar. In fact there was a pro- 
gressive diminution in the value of E;—E; when the concentration of iron 
(III) was increased from 0-3 to 1-5 millimolar. 


TABLE III 
Effect of sodium carbonate on the polarographic behaviour of iron (IIT) 
in presence of triethanolamine 








Conc. of Conc. of Conc.of  (E,);* (Ej). * 

Triethanol- Sodium Iron VS. vs. (Ez; —Ey) (E; —Epe 
amine carbonate (Milli- S.C.E. S.C.E. 

(Molar) (Molar) molar) (Vol.) (Volt) (Volt) (Volt) 
0-02 1-00 1-182 —0-843 —1-531 0-061 0-076 
0-02 1-00 3-546 —0-799 —1-540 0-064 0-074 
0-10 1-00 1-182 —0-867 —1-557 0-066 0-077 
0-10 1-00 5-910 —0-813 —1-546 0-068 0-079 
0-10 0-10 1-182 —0-774 —1-504 0-072 0-075 
0-10 0-10 5-910 —0-760 —1-553 0-070 0-029 
1-00 1-00 1-182 —0-906 —1-637 0-067 0-070 
1-00 1-00 5-910 —0-874 —1-637 0-072 0-075 
1-00 0-10 1-182 —0-820 —1-635 0-078 0-068 
1-00 0-10 5-910 —0-785 —1-649 0-092 0-100 





* Ey represents the half-wave potential and the subscripts 1 and 2 refer to the first and 
the second waves respectively. 


In 0-02 M triethanolamine with an increase in the concentration of 
iron, the first wave is divided into two poorly separated waves, the second 
one being much smaller than the first (Fig. 3). The complexity of the wave, 
however, does not affect the constancy of the value of the diffusion current 
constant. 


(c) Ammonium chloride (1 M)—ammonium hydroxide (1 M).—Some of 
the polarograms showing the effect of concentration of triethanolamine 
are given in Fig. 4, The results are summarized in Table IV. It is clear 
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that the half-wave potentials progressively increase with increasing 
concentration of triethanolamine and finally attain a constant value. The 
first wave is reversible while the second wave is irreversible. It has been 
noticed that the height of the second wave is double the first in all the cases, 


TABLE IV 


Effect of ammonium chloride (1 M)—ammonium hydroxide (1 M) on 
polarographic waves of iron (III) in presence of triethanolamine 


Conc. of iron 1-182 millimolar 








Conc. of Ist wave 2nd wave 
triethanol- Half-wave (E;—E}); Half-wave (E;—E,), 
amine potential potential 
(E,), vs. S.C.E. E(y)2 vs. S.C.E. 

(Molar) (Volt) (Volt) (Volt) (Volt) 
0-02 —0-449 0-058 —1-436 0-051 
0-10 —0-491 0-056 —1-463 0-043 
0-30 —0-519 0-062 —1-468 0-043 
0-50 —0-528 0-057 —1-476 0-040 
0-70 —0-535 0-062 —1-483 0-043 
1-00 —1-532 0-056 —1-480 0-038 





3. Effect of concentration of iron (III) on the diffusion current in base solutions 
containing triethanolamine 


Table V shows the effect of concentration of iron on the diffusion current 
constant under varying conditions. The diffusion current constant only 
for the first wave has been given, since the second wave is not suitable for 
analytical purposes. The results indicate that in base solutions containing 
sodium carbonate the diffusion current constant decreases with an increase 
in the concentration of iron except in solutions containing 0-02 M triethanol- 
amine + 1M sodium carbonate. In these solutions the solubility of iron 
is limited. Hence the presence of sodium carbonate is not of any advantage 
in the polarographic estimation of iron. The most suitable base solutions 
for the estimation of iron contain 0-1 M triethanolamine + 0-1 to 1-0 
molar sodium hydroxide but in presence of ammonium chloride-ammonium 
hydroxide, the concentration of triethanolamine has to be 1-0 molar due to 
the low solubility of the complex. 
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TABLE V 
Effect of concentration of iron (III) on the diffusion current constant 


The values mt ttin mg.? sec.~+ has been given at —1-4V, —1-1 V and—0-85 V vs. S.C.E., 
in solutions containing sodium hydroxide, sodium carbonate and ammonium chloride-ammonium 
hydroxide respectively. In this table B.S., X, and Trien represent base solution, mt ft, and 


triethanolamine respectively. 








Conc. of Diffusion 
ferric iron current 

(Millimolar) constant 
lg 

Cmt tt 


Conc. of Diffusion 
ferric iron current 
(Millimolar) constant 
_ 
Cm tt 





(1) B.S.: 0-02M Trien-+-0-1M NaOH 


X = 1-361 
0-296 "1-95 
0-591 1-95 
1-182 1-85 
2-364 1-84 
4-728 1-84 
(2) B.S.: 0-1M Trien + 0-1M NaOH 
X = 1-361 
0-296 1-86 
0-591 1-85 
2-360 1-87 
4-720 1-86 
7-080 1-85 
(3) B.S.: 1-00M Trien + 0-10M NaOH 
X = 1-379 
0-296 2-14 
1-184 1-70 
2-368 1-69 
8-274 . 1-57 
16-548 1-51 
(4) B.S.:0-02M Trien + 1-00M NaOH 
X = 1-589 
0-380 1-75 
1-520 1-72 
7-980 1-71 
(5) B.S.:0-10M Trien + 100M NaOH 
X = 1-382 
0-296 1-94 
1-184 1-71 
5-920 1-69 


11-840 1-70 





(6) B.S.:0-002M Trien + 1-00M NaOH 


X = 1-580 
0-448 1-07 
1-120 1-07 
2-240 1-04 
(7) B.S.:0-02M Trien + 5-:00M NaOH 
X = 1-579 
2-240 1-00 
5-600 0-99 
(8) B.S.:0-10M Trien + 5-00M NaOH 
X = 1-577 
0-448 0-99 
1-120 0-99 
2-240 0-98 
5-600 0-97 
(9) B.S.:0:02M Trien +1-00M Na,CO, 
X = 1-415 
0-473 1-23 
2-364 1-21 
(10) B.S.:0-10M Trien+1-00M Na,CO, 
X = 1-391 
1-182 1-35 
5-910 1-27 
(11) B.S.:1-00M Trien+-1-:00M Na,CO, 
X = 1-379 
1-182 1-30 
5-910 1-17 
(12) B.S.:0-10M Trien+-0-10M Na,CO, 
X = 1-607 
1-158 1-60 
5-790 1-33 
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TABLE V—({Contd.) 





Conc. of Diffusion Conc. of Diffusion 

ferric iron current ferric iron current 

(Millimolar) constant (Millimolar) constant 
_ig ig 

Cmi tt Cmi tt 


(13) B.S.:1-00M Trien+0-01M Na,CO, (14) B.S.: 1-0M Trien + 1M NH,Cl 


X = 1-602 + 1M NH,OH 
1-158 1-44 . ¥=1-427 
2-316 1-36 0-306 1-86 
5-790 1-26 1-222 1-84 
11-580 1-19 2-444 1-87 
6-110 1-82 














The use of 5M sodium hydroxide decreases the diffusion current con- 
stant and the solution becomes viscous without any obvious advantage. It 
may also be pointed out that in mosi cases the value of the diffusion current 
constant is slightly higher at lower concentrations of iron (below 0-3 milli- 
molar). 


4. Use of sodium acetate 





In these experiments 0:2 M sodium ions have been included to eliminate 
the migration current. The second wave is too very close to the final current 
rise. At pH 8-3 the half-wave potentials for the first wave of iron are 
— 0-413 V, — 0-422 V, — 0-440 V, and — 0-453 V, vs. S.C.E., when the base 
solutions contain 0-02, 0-10, 0-50, and 1-00 molar triethanolamine res- 
pectively. The values of E; — E; are 0-080V, 0-086 V, 0-060 V, and 
0-060 V indicating that the wave is reversible at higher concentrations of 
triethanolamine. Due to solubility difficulties the effect of concentration 
of iron on the diffusion current has been tried employing 0-5 M triethanol- 
amine solution. These solutions have an yellow colour at higher concentra- 
tions of iron, perhaps due to the formation of colloidal hydrous ferric oxide. 
The diffusion current constant decreases by about 18% when the concentra- 
tion of iron is changed from | to 5 millimolar at pH 8-3 and also at pH 7. 











(B) Base Solutions Containing Diethanolamine 


The behaviour of diethanolamine complexes of iron resemble closely 
the triethanolamine complexes but the solubility of the complex is very much 
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less in diethanolamine solutions. Table VI gives the results obtained with 
sodium hydroxide and sodium carbonate. The other base solutions employed 
in the case of triethanolamine have not been used due to the low solubility 
of the complex. The half-wave potential of the first wave is very nearly 
constant at different concentrations of the amine while that of the second 
wave increases with an increase in the concentration of the amine. The 
E; — E; value of the second wave varies with the variation in the concen- 
tration of the amine reaching almost a constant value at higher concentra- 
tions of the amine. In these base solutions the solubility of the diethanol- 
amine complex is quite low for the method to be of practical value in esti- 
mating iron. 
TABLE VI 


Effect of sodium hydroxide and sodium carbonate on the polarographic waves 
of iron (III) in presence of diethanolamine 





Composition of (Eye 
the base S. (E vs. (E; —E,) 
solution EB. S.C.E. 
(Volt) (Volt) (Volt) (Volt) 





(1) 0-10M NaOH + 0:30M Dien... —1-026 0-075 —1-576 0-064 
(2) ys +0:65M ,, .. —0-938 0-078 —1-611 0-076 
(3) +1: ‘00M ,, .. —0:928 0-069 —1-624 0-076 
(4) 1-00M ‘NaOH +0:30M ., .. —1:005 0-075 —1-637 0-048 
(5) es 4 0:65M .. —0:995 0-065 —1-675 0:055 
(6) - * 100M ,, .. —0:994 0-062 —1-709 0-065 
(7) 500M NaOH + 0:30M ,, .. —1:047 0-069 —1-639 0-035 
(8) oe + 1:00M .. —1:040 0-059 —1-681 0-038 
(9) 0-10 M Na,CO, + 0-30 Dien .. —0-821 0-086 —1-466 0-044 
+ 0:005% gelatin 
(10) - + 0:65M Dien —0-834 0-064 -493 0-051 
+ 0:005% gela- 
tin 
(11) - + 1-00M Dien —0-830 0-064 “512 0-049 
+ 0:005% gelatin 
(12) 1-00M Na,CO, +0-30M Dien —0°816 0-074 -519 = 0-062 
+0-005% gelatin 
(13) ji +0-65M Dien —0-833 0-068 ‘535 = 0-048 
+0-005% gelatin 
(14) m +1-00M Dien —0-837 0-060 ‘550 0-048 
+0:005% gelatin 





It has also been noticed in the present investigation that at high con- 
centrations, the ratio of the height of the second to the first wave decreases 
with increasing concentrations of the amine. The effect, however, is not 
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sO prominent as in the case of triethanolamine solutions. In base solutions 
containing 0-10 M sodium hydroxide the ratio changes from 1-9 to 1-6 
when the concentration of diethanolamine is changed from 0-3 to 1-0M. 


(C) Base Solutions Containing Monoethanolamine 


The solubility of the monoethanolamine complexes of iron is very low 


in all the base solutions that have been used in the present work. Hence 
no detailed studies have been made. 


DISCUSSION 
(A) Base solutions Containing Triethanolamine 
1. Effect of pH on the half-wave potential 


Fig. 5 indicates the variation of half-wave potential with pH. The 
values of E; — E; given in Table I indicate that the ferric to ferrous reduction 
is irreversible between pH 11 and 8. The dotted line indicates the reversible 
half-wave potential calculated by employing Eyring’s method® and it can be 
seen that the form of the curve remains the same. It can therefore be con- 
cluded that the different slopes are due to changes in the nature of the com- 
plexes that are being reduced at the dropping mercury electrode. The 
region AB corresponds to the complex of triethanolammonium ion with 
ferric iron. The low value of the slope at BC is perhaps due to a compa- 
ratively slow change in the ratio ion/amine around the pK value of the 
triethanolamine. The high value of the slope of the curve CD indicates 
that ferric-triethanolamine complex is getting formed. The change of the 
slope at DE is indicative of partial modification of the complex by the entry 
of hydroxyl group. Méeites? has also suggested that the hydroxyl group 
gets into the complex at high pH values. When the strength of the alkali 
is increased considerably, the introduction of the hydroxyl group into the 
molecule probably ceases and the hydroxo-triethanalo complex remains 


unaltered with the result that in the region EF there is very little change in 
the half-wave potential. 


2. Effect of pH on the ratio of heights of the two waves 


It has been noticed (Tables I and I) that at high concentrations of the 
alkali, the ratio of the heights of the second to the first wave is not equal 
to two but varies from nearly 2 to 0-2. Since the second wave is a two 
electron process and the first a one electron process it is reasonable to expect 
(according to Ilkovic Equation®) the ratio to be equal to two. Similar low 
values for the ratio have also been obtained by Komarek® in alkaline 
mannitol solutions of iron (III). The above remarkable fact can be explained 
by assuming the existence of two types of ferrous complexes, the half-wave 
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potential of one being more positive than the other and the equilibrium 
between the two complexes being slow at the dropping mercury electrode. 
The first wave is due to the complex getting discharged at a less negative 
potential. The wave having a more negative potential does not appear 
as it gets masked with the final rise of current due to sodium discharge. 
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3. Nature of complex formation 


In the case of reversible reduction to the metallic state, the half-wave 
potential of the complex metal ion is related to the excess concentration 
of the complexing ion according to the equation’® 
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AF, 


A log Cy waktiel (0-06014/n) at 30° Cc. (1) 


where q is the number of ligands attached to the complex metal ion, n the 
number of electrons involved in the reduction process, Cy the excess con- 
centration of the complexing agent, and 4 E, the shift in the half-wave 
potential. 


In the case of reduction from one oxidation state to another, the above 
equation can be modified as 


Floats ~~ @ — 9 0-06014/a) 2) 


where p is the number of ligands attached to the oxidised complex, q the 
number of ligands attached to the reduced complex, and a the number of 
electrons involved in the reduction process. 


For irreversible reductions the following equations can be written based 
on the treatment given by Tamanushi and Tanaka." 


4E; 


Flog C, = ~ 4 (0-06014/an) (3) 


where a = the fraction of the total applied potential that favours the for- 
ward reaction. 


The value of a is obtained from the following expression :— 
E; — E; = — (0-0574/na) (4) 


In the case of reduction from one oxidation state to another, equation 2 
can be modified as 


ae abe = — p(0-06014/aa) (5) 

It is interesting to note that in irreversible reductions from one oxida- 
tion state to another the number of ligands attached to the reduction com- 
plex does not enter the equation. The application of equations 1 and 3 
for the reduction of ferrous ion to metal seems to be justifiable since freshly 
deposited iron is in the form of amalgam.'? 


The application of the above equations to obtain the formule of the 
complexes formed is beset with many difficulties: (1) the waves show vary- 
ing degrees of irreversibility with a change in the concentration of the com- 
plexing agent and the equations for the irreversible wave can be applied only 
when the value for the irreversibility is constant; (2) the ratio of the second 
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to the first wave of iron complexes in sodium hydroxide solutions is much 
less than two; (3) in sodium carbonate solutions, when 1-00 M triethanol- 
amine is present, the half-wave potential of the second wave does not change 
when the concentration of sodium carbonate is changed from 0-10 to 
1:00M; a definite shift however, is noticed when 0-1 M triethanolamine 
is present; the half-wave potential of the first wave is susceptible to changes 
in the concentration of iron; further, the irreversibility of the first wave is 
not sufficiently high to warrant the elimination of g in p — q of the reversible 
wave equation. 


(a) Base solutions containing sodium hydroxide and triethanolamine.—Ref. 
to Table II indicates that (E;), is practically constant at low concentrations 
of triethanolamine and the wave is reversible.’° Hence it has to be concluded 
that p = q under these conditions. Equality of p and q is also established 
for all concentrations of triethanolamine when the base solution has large 
amounts (5M) of sodium hydroxide. With higher concentrations of tri- 
ethanolamine, the half-wave potential of the complex diminishes. This is 
rather surprising since it is usual to find an increase in the value of the half- 
wave potential. Similar behaviour has also been noticed by Meites™ in 
systems containing copper-ammonia-citrate complex but a suitable explana- 
tion is not forthcoming at present to explain such exceptions. 


(b) Base solutions containing sodium carbonate and triethanolamine.— 
At a given pH the half-wave potentials in sodium carbonate base solutions 
have been found to be 0-1 V more positive than in base solutions contain- 
ing sodium acetate, indicating the formation of new type of complex (the 
carbonato-ethanalo complex). It is difficult to get the formula of ‘the 
ferric complex. Applying equations 3 and 4 the formula of the ferrous 
complex can be written as follows :— 


Fe (II) (Trien), (CO,);. 


(c) Base solutions containing ammonia-ammonium chloride and tri- 
ethanolamine.—In this base solution the first wave is reversible while the 
second wave is slightly irreversible. At high concentrations of triethanol- 
amine (above 0-7 molar) the half-wave potential of the first wave is almost 
a constant indicating that p=q. Applying equations 2, 3 and 4 for the 
concentration range of triethanolamine between 0-1 and 0-7 molar one 
gets p—q=1 and q=1. The following formule can be written for the 
complexes : 


[Fe (III) (T rien), (NHs) }***; [Fe (II) (Trien); (NHs)]}* 
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(B) Base Solution Containing Diethanolamine and Sodium Carbonate 


In this base solution the first wave is reversible at higher concentra- 
tions of diethanolamine. The second wave, however, is irreversible and 
the irreversibility reaches a constant value at higher concentrations of the 
amine. Applying equations 2, 3 and 4 the following formule can be written 
for the complexes formed: 


[Fe (III) (Dien), (CO;),]* and [Fe (II) (Dien), (CO,),]. 


SUMMARY 


1. Employing ethanolamines as the complexing agents a systematic 
study on the polarographic behaviour iron has been carried out at various pH 
values and in the presence of the following electrolytes: (a) sodium hydroxide, 
(b) sodium carbonate, (c) acetate, and (d) ammonium chloride-ammonium 
hydroxide. 


2. In triethanolamine solutions containing sodium hydroxide, sodium 
carbonate, and ammonium chloride-ammonium hydroxide, the polaro- 
grams consist of two waves corresponding to the reduction of ferric to ferrous 
and ferrous to iron. The ratio of the heights between the second and the 
first wave is less than two in the case of sodium hydroxide while in the other 
two cases it equal to two. The variation of half-wave potential with pH 
has indicated the formation of complexes involving triethanolammonium 
ion, triethanolamine, and hydroxyl group. 


3. The suitability of using triethanolamine in base solutions for polaro- 
graphic estimation of iron has been discussed. 


4. Experiments conducted with mono- and diethanolamines show that 
the solubility of the complexes is very much lower in the base solutions which 
makes the study very difficult. 


5. The formation of the following complexes have been indicated from 
the polarographic measurements: 


(1) [Fe (II) (Dien), (CO;),]* (2) [Fe (II) (Dien), (CO,),] 


(3) [Fe (11) (Trien), (CO;);] (4) [Fe (IIT) (Trien), (NH;)z}*** 
(5) [Fe (11) (Trien); (NH5)y}** 
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Effect of sodium hydroxide (0-1M) onthe polarographic behaviour of iron (1+182 millimolar) 
in presence of triethanolamine 
Curves 1, 2, 3 and 4 have been taken using 0-002, 0-02, 0-10 and 1-00 M triethanolamine. 


Sensitivity: s/30, s/50, s/50, s/50 respectively. 
The voltage line just before the beginning of the polarogram corresponds to —0+573 V ws. S.C.E. 
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FIG. 3 
Effect of sodium carbonate (1 M) and triethanolamine (0-02 M) on the polarographic behaviour 
of iron. 
Curves 1 and 2 have been obtained with 1-182 and 3:546 millimolar iron at s/50 and s/150. 
Curves 3 and 4 are the blanks at s/50 and s/150. 
The voltage line just before the beginning of the polarogram corresponds to —0+473 V vs. S.C.E 
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Effect of ammonium chloride (1 M) and ammonium hydroxile (1 M) on the polarographic behaviour 
of iron (1-182 millimol.r) in presence of triethanolamine. 

Curves 1, 2, 3 and 4 have been taken with 0°02, 0°10, 0°30 and 9-50M_ triethanolamine at s/30, 
s/50, s/50 and s/50 respectively. 


There has been slight precipitation under the conditions at which 
curve | is recorded. 


The voltage line just before the beginning of the polarogram corresponds to —0+273 V ws, S.C.E, 
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